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CXXXVIIT. The Thickness of the Saturated Helium II Film’ 


By Raymonp Bowers 
Institute for the Study of Metals, University of Chicago* 


[Received August 18, 1953] 


ABSTRACT 


The average thickness of the saturated He IT film on an aluminium foil 
has been measured by a gravimetric method. For the film in contact with 
the bulk liquid, values between 120 and 180 atomic layers (4:8 x 10-* cm 
and 7-2 10~-* cm) have been observed, depending on the height of the 
specimen above the bulk liquid level. The average thickness is indepen- 
dent of temperature except within 0-006 degrees of the lambda point. 
The values obtained suggest a film profile d=11-8x10~*/h1/2, where d 
is the actual film thickness in cm at a height h above the liquid level. 
For the He II film not in direct contact with the liquid, a thickness of 
55 atomic layers (2-2 x 10~® cm) has been observed, independent of temper- 
ature and position of the bulk liquid level. This value is discussed in terms 
of a lack of temperature equilibrium. Observations of the evaporation 
of the film by radiation influx have been made. The possibility that the 
disparity between the results of various film determinations may be due to 
differing conditions of heat influx is considered in the discussion. 


§ 1. INTRODUCTION 


In spite of numerous measurements of the thickness of the saturated 
helium II film, this fundamental property of the film is still not known with 
any certainty. The earliest experiments of Daunt and Mendelssohn (1939) 
and those of Kikoin and Lasarew (1938, 1939) indicated that the saturated 
film was approximately 100 atomic layers. Actually, the individual results 
at different temperatures obtained by Daunt and Mendelssohn varied in a 
random manner from 2:7 to 4:910-§cem (70 to 120 atomic layers). 
More recent work has yielded a considerable amount of data for the thick- 
ness of the film, but there is little consistency between the data of different 
workers. In table 1 we have summarized the principal results of a number 
ofinvestigations. For simplicity, we shall consider the general magnitudes 
Pip aay Eds Ae aod it AO SS GO AES AE TER Se eS a ee 
ye i uthor. 

tT eat ie or a aaa ae in the paper refers to the film in direct contact 
with the liquid.. The special case of a film isolated from the liquid except via 
the vapour is discussed in the text, but it is referred to specifically as the 
‘isolated film ’. 

4Q2 
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observed ; and we have not included fine details such as variation with 
temperature and height of the film above the bulk liquid level. A more 
detailed analysis of previous work is given by Atkins (1950). 


Table 1 
Investigation Method Substrate Thickness 
Daunt and Mendelssohn 1939 Volumetric Metal 3-6 x 10-* cm 
Kikoin and Lasarew Volumetric Metal 2-3 x 10-® em 
1938, 1939 
Burge and Jackson 1949 Optical Thin layer 2x10-* cm, varying 
of barium with height and 
stearate slightly with tem- 
perature 
Long and Meyer 1949 Volumetric Rouge 6-4 x 10-® cm 
absorption 

Atkins 1950 Film oscilla- Glass 2x 10-% cm, varying 
tion with height, tem- 
perature and 


nature of surface 


(Helium film of 100 atomic layers is equivalent to 4 10-® cm) 


It is evident from the table that there are large disparities between the 
results of various investigations. Actually, they are larger than is 
immediately obvious from the table. For example, the result of Burge and 
Jackson is the thickness at a height of 1 cm above the liquid level ; on the 
other hand, Daunt and Mendelssohn’s value is an average over 3 to 18 cm 
of film above the liquid level. If the thickness (d) of the film varies with 
height (2) above the bulk liquid level according to the relation 


d=k/him er Eh ay 


where k is a constant and n=2 or 3 (Frenkel 1940, Schiff 1941, Bijl, de 
Boer and Michels 1941) then an extrapolation of the Burge—Jackson result 
to the heights involved in the experiment of Daunt and Mendelssohn 
results in a discrepancy of about five times. 

It is not obvious that there are systematic differences between the 
various investigations, of a type which might explain the gross disparities. 
Tn allied studies, such as those of adsorbed films, one would immediately 
consider the nature of the surfaces used. However, apart from consider- 
ations such as those of Frenkel (1940) involving surface energies and those 
of Bowers (1953 b) for the special case of a fine powder, it is not easy to see 
how the substrate could determine the thickness of a film as thick as the 
one being considered here. The forces between a helium atom and the 
substrate must be of very short range compared with any of the measured 
values for the Helium II film thickness (see case of classical liquids, 
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Bowden and Throssell 1951). The extent to which a substrate can affect 
the film is best seen from the value of H,, the energy of adsorption of the 
first atomic layer. The influence of substrate on E, has not been studied 
extensively for helium and it is complicated by the fact that H, varies 
with amount adsorbed (Keesom and Schweers 1941), but it is well known 
that for other gases HZ, varies very little as the substrate (non-ionic) is 
changed (Brunauer 1945, p. 244). 

In considering a possible variation of film thickness with substrate, it is 
tempting to use the observed independence of film transfer on the chemical 
nature of the substrate (Chandrasekhar and Mendelssohn 1952) as an 
argument for a film thickness which does not change with surface. 
However, the transport property is not a good criterion for considering the 
thickness of the film. The total mass transfer is the product of a linear 
velocity in the film and the thickness involved and there is some reason 
for believing that this product is invariant with respect to a change in 
linear velocity or thickness (see Dingle 1952, Equation 135). 

The discussion given above indicates some of the difficulties in under- 
standing the disparity between the results of different investigations. The 
investigation to be described in this paper was undertaken to obtain 
further data for a resolution of these difficulties. The film thickness has 
been measured by a gravimetric method. This method is direct, there 
being no corrections necessary to the observed weight of the helium IT 
film apart from the evaluation of the surface roughness factor of the 
specimen. The technique used was the same as for similar studies of 
unsaturated films (Bowers 1953 a, 1953 b). The microbalance assembly 
has been described elsewhere (Bowers and Long 1953). A summary of 
the results has already been published (Bowers 1953 c). 


§ 2. MrerHop 


The measurements were made on an aluminium foil 7-6 cm x 5 cm and 
0-002 cm in thickness. The shape of the specimen is of great importance 
in this investigation since it must have no folds or pores of a type that would 
retain bulk liquid. Since the thickness of the film forming on it is deter- 
mined from the effective increase in weight of the specimen, it is essential 
that no drops of bulk liquid be included in the total weight. The form 
shown in fig. 1 was the one found most suitable for our work. The foil 
was folded into a series of smooth ‘s’ bends. Care was taken to ensure 
that the folds of the foil did not touch each other and that no sharp bends 
appeared in it. 

In order to support the specimen from the wires which hang from the 
beam of the microbalance, holes were cut into the top of the foil. The 
holes had a diameter of 1 mm and were large enough so that a sheet of 
liquid would not be maintained across them by surface tension. A narrow 
tail (diameter 0-002 in.) was attached to the bottom of the specimen in 
order that contact could be made between the specimen and the liquid 
below. 
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The design of a suitable specimen has been the greatest single difficulty 
in this investigation. For our work on unsaturated films, it was sufficient 
to use a loosely wound coil, but at saturation the formation and retention 
of liquid between the folds rendered a precise determination of the film 
thickness impossible. Experiments at saturation, using such a specimen, 
yielded very high values for the film thickness with little consistency both 


Fig. 1 


or 
ALUMINUM 
SPECIMEN poem 
Ho 
WIRE TAIL 2.7¢m 
(DIAM. .002" ) 


Hy 


above and below the lambda point. It was observed that the apparent 
thickness could be increased by tightening the windings of the coil, 
indicating the presence of bulk liquid. After numerous experiments with 
specimens of this kind, it was concluded that the specimen must approxi- 
mate more closely the conditions of a plane sheet. In using the specimen 
shown in fig. 1, it was necessary to reduce the area of the foil by a factor 
4 since it is not as compact as the coils mentioned. Hence a fourfold 
decrease in the sensitivity was expected. 

The specimen was cleaned and outgassed as in the work on unsaturated 
films. Its effective area was determined from a measurement of the 
nitrogen isotherm. The effective area was found to be 97 cm?, giving a 
roughness factor of 1-3, in agreement with our previous work. Thus 1 
atomic layer of helium on the specimen weighs 0-40 g, assuming the 
liquid spacing of 4.4. (The balance sensitivity is 0-0158 mm/ug.) 

With the exception of the specimen and some alterations in the radiation 
shielding, the arrangement used was precisely the same as that used in the 
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earlier work (Bowers 1953 a, 1953 b). We have improved the radiation 
shielding by placing an extra shield on the specimen support wire itself at 
10 cm above the specimen. (See fig. 4, Bowers and Long 1953.) 

The liquid in the specimen chamber is condensed in at a slight over- 
pressure. While this is being done, drops form on the specimen, but they 
appear to evaporate or fall off in a few minutes ; the fact that measurements 
of the saturated film of He I and also that of liquid nitrogen gave about 
10 atomic layers confirms our belief that no appreciable amount of bulk 
liquid adhered to the specimen. In the case of He II, the measurements 
made were as follows: the change in weight of the specimen was observed 
as the liquid helium level (H) in the specimen chamber was raised from a 
height H, at the bottom of the chamber to a height H, within 3 mm of the 
bottom of the foil. It can be seen from fig. 1 that until the liquid level 
touches the tail of the specimen, the liquid and the foil are not in direct 
contact ; on touching the tail, contact is made via the mobile film on this 
tail. The zero position for the balance is determined as in our previous 
work by observing the position of the balance beam at 15°% saturation. 

The film thickness obtained on dividing the change in weight of the 
specimen by the product of its area and the density of the film, is an 
average thickness over the specimen. In considering the values obtained, 
we must take account of the film profile, i.e., the variation of the actual 
film thickness with height above the liquid level. If the profile of the 
film is given by eqn. (1) above, then the weight of film is obtained by 
integrating (1) over the height of the specimen. It is easily shown that 
the weight w is given by 


eer Pm ee, 2... (2) 


where p is the density of the film, A is the total effective area of the foil 
(97 cm?) and a is the vertical length of the foil (7-6cm). The limiting 
heights 4, and h, are shown in the inset of fig. 2. Expressed in another 
form, the average number of layers (1) 


Palas 
pAt n—1 at 
where t is the atomic spacing (4 A for all but the first few layers). Hence 
the film profile can be decided from our observations. 


Ni i : ; . i (3) 


§ 3. RESULTS 


The first measurements made with liquid helium were to check our 
data on the unsaturated film. The principal observations were repeated 
and verified ; close to saturation (p/py=0-9992) the film is only 20 atoms 
thick. When this had been established, liquid helium II was condensed 
into the specimen chamber and the average film thickness was observed. 
as the liquid helium level was raised from the bottom of the specimen 
chamber (H, in fig. 1) until it touched the tail of the specimen. This 
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thickness of the isolated film was found to be 55 atomic layers, independent 
of temperature and also independent of the level of the liquid. 

These observations will be discussed later, but it is convenient to men- 
tion some aspects of the observations which concern the sensitivity of our 
method. In our work on the unsaturated film, using a specimen of four 
times the area, it was possible to measure thicknesses to within one atomic 
layer. In this case, the accuracy seemed to be limited by the ultimate 
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of the bulk liquid level. 


sensitivity of the balance. However, when measuring the film thickness 
in the presence of bulk liquid, we have found the scatter was higher than 
the four atomic layers expected. When liquid is present, the balance is 
less stable ; small fluctuations in temperature produce large spurious 
deflections of the balance. Even when controlling the temperature to 
1/1000 of a degree, we have not found it possible to reduce the scatter 
below 5-10 atomic layers. The greater instability observed when bulk 
liquid is present is probably due to the fact that when a slight temperature 
change occurs, a small amount of helium evaporates or condenses in the 
liquid below the specimen. The resulting movement of gas disturbs the 
balance. This type of instability might ultimately limit the extent to 
which the sensitivity of the gravimetric method can be increased. 

We shall now describe our observations of average film thickness as the 
liquid level is raised further so that it is in contact with the wire tail. 
The variation of this thickness with temperature for various positions of 
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the bulk liquid level in contact with the tail is shown in fig. 2, (Small 
corrections have been made for the effect of buoyancy and surface tension 
on the immersed part of the tail.) In this figure, h, is the distance of the 
bulk liquid level from the bottom of the foil. It can be seen that the values 
illustrated refer to distances of the liquid level from the bottom of the foil 
which vary between 2-6 cm and 0-3cm. The experimental points have 
been connected by straight lines to distinguish corresponding data. 
We believe that the thickness, for a given height of the bulk liquid, is 
independent of temperature within our experimental error. 


Fig. 3 
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With regard to the film profile, our data are best described by taking 
n=2 in eqn. (3) for all temperatures. In this case 
25 
N= A eee ee Rc ge Fos ~ yt, (4) 
Graphs of N against h,/2—h,/? are shown in fig. 3 for a number of 
temperatures. The values of K calculated for different temperatures are 
shown in fig. 4. This graph, which is really a smoothed representation of 
all the data, shows the lack of dependence of film thickness on temperature 
(except very close to the lambda point). 
Observations have also been made on the thickness of the saturated film 
above the lambda point. It is more difficult to measure the film thickness 
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above the lambda point because the temperature cannot be held very 
steady. At 2-25°K we estimate the thickness of the saturated film to be 
between 10 and 15 atomic layers. This agrees with the work of Jackson 
and Henshaw (1953), but is smaller than the values due to Kistemaker 
(1947). We have also attempted to observe the precise temperature 
interval in which the thickness of the He IT film falls from approximately 
150 layers to the 10 layers observed above the lambda point. For this 
observation our specimen has certain disadvantages compared with that 
used by Burge and Jackson. The rate at which the film forms on our 
specimen is governed to some extent by the film creeping up the wire tail 
(diameter 0-002 in.). At the lowest temperatures, it is easily calculated 
that creep up the wire alone could furnish the complete film in about ten 
minutes. (Of course, equilibrium is reached in a shorter time.) However, 


Fig. 4 
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very close to the lambda point, film creep up the wire becomes very small 
and above 2:17°K it is necessary to wait for more than a half hour for 
equilibrium. It is not easy to increase the diameter of the tail because 
buoyancy and surface tension corrections for the immersed part become 
large. . 

Taking care to ensure that equilibrium had been reached, we found that 
the film maintained its full thickness to within 0-006 degrees of the lambda 
point. It was difficult to place more precise limits on the temperature 
interval than this because the experiments become very tedious nearer 
the lambda point. Our observation is essentially in agreement with that 
of Burge and Jackson who quote the interval as 0-004°K. 

A series of observations have been made with the purpose of evaluating 
the influence of heat influx by radiation on the film thickness. It will be 
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more convenient to discuss the results recorded so far before describing 
these observations. 


§ 4. Discussion 


We shall first consider the case of the film which is isolated from the 
bulk liquid. This film was observed to be thinner than the film in contact 
with the liquid by a factor between 2-5 and 3-3. The heat contact due to 
the film creeping up the tail is not present in the case of the isolated 
film ; accordingly, we interpret the reduction in film thickness as being 
due to a lack of complete temperature equilibrium. An experiment by 
Daunt and. Mendelssohn (Mendelssohn 1946) gave the result that an 
isolated film under conditions similar to these was only one-tenth as thick 
as when the specimen was in contact with the bulk liquid. The difference 
between their 10% and our 30% is attributed to fine differences in heat 
influx and the actual values obtained have little significance since they 
refer to a steady state but not to temperature equilibrium. Our inter- 
pretation of the smaller thickness of the isolated film is supported by the 
observation that for the film in contact with the liquid, the approach to 
equilibrium is slow very near the lambda point where the film creep is low. 
Nevertheless, it is surprising that heat and mass exchange between the film 
and the ambient vapour should be small enough for the creep rate to have 
this influence on equilibrium time. It is easily calculated that under the 
equilibrium conditions of our experiment, about 1-510? molecules 
(1 g) hit 1 cm? of foil per second (Mayer and Mayer 1940). With such a 
large flux, only a small accommodation coefficient between vapour 
molecules and the film can explain our observations.* 

The thickness of the film in contact with the liquid measured in this 
investigation is more in accord with the results of Daunt and Mendelssohn 
and also of Long and Meyer than with those of Atkins and of Burge and 
Jackson. Actually, if our result d=11-8 x 10-®/h1/? is averaged over the 
heights involved in the Daunt and Mendelssohn experiment, close 
agreement is obtained with their data. The value of the constant K 
observed by us is about 6 times that given by Atkins and by Burge and 
Jackson. However, the value of K recorded by Burge and Jackson arises 
from measurements of film thickness at about 1 cm above the liquid. It 
would not be correct to assume that our value of K gives the actual 
thickness of the film at a height of 1 cm above the liquid, because this 
would involve an extrapolation of our data well into a region which is not 
measured in the experiments. We have not checked the validity of 
eqn. (1) at small heights ; the highest film thickness observed by us refers 
to a film extending from 0-3 em to 7-9 cm above the liquid level. 

In looking for a general scheme into which one could fit the variety of 
values recorded by different workers, it is well to observe that, with the 
exception of the value due to Long and Meyer, the film thickness appears 


* We are indebted to Professor L. Meyer for drawing our attention to this 
point. : : 
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to be greater on metal surfaces than it is on dielectrics, as suggested by 
Frenkel. The data of Long and Meyer were obtained from an adsorption 
measurement on finely divided rouge. It has already been pointed out 
(Bowers 1953 b) that one would expect very high values for the film 
thickness in this case, since the measurement corresponds to a film at an 
effective height of the order of the grain size and this is very small. 
However, as a method of measuring the saturated film thickness, as opposed 
to the adsorption isotherm, the method of Long and Meter is subject to 
doubt. It is difficult to be sure that capillary condensation was not 
occurring very close to saturation or even some hysteresis due to over- 
saturation in the powder. Of course there is still insufficient data avail- 
able to decide whether a scheme based on Frenkel’s suggestion is valid. 
Indeed, reasons were given in the introduction to this paper why one 
would not expect the chemical nature of substrate to affect the film 
thickness. 

In comparing the results of previous investigations, we are hampered 
by the fact that no values are given for the surface roughness factor of the 
specimens used. It is well known that for some metal foils the effective 
surface area for adsorption may be many times the apparent geometrical 
area (Brunauer 1945, p. 284). This is important in evaluating the results 
of Daunt and Mendelssohn and also of Kikoin and Lasarew, although the 
direct optical determination of film thickness by Burge and Jackson does 
avoid this difficulty. A high roughness factor in the specimens used by 
the former could provide an explanation of why their recorded values for 
the film thickness are higher than those of Burge and Jackson. Such an 
explanation would not, however, cover the values recorded in the present 
investigation, in which the roughness factor has been taken into account. 

The value of m in eqn. (1) was found to be 2 for our data. This is the 
value predicted by Bijl, de Boer, and Michels (1941) on the basis of a theory 
involving zero point energy. A value n=3 as suggested by Frenkel 
(1940) and by Schiff (1941) cannot be fitted to our data. This is of 
interest because the calculation of Bijl, de Boer, and Michels is not 
satisfactory since it involves unrealistic forms of the wave functions in 
the film (Mott 1949). 

The calculations of Frenkel and of Schiff were essentially a consideration 
of the film which would result from the van der Waals attraction to the 
wall. The profile given by them resulted from the effect of gravity on the 
film system. They suggested that the essential difference between He II 
and He I, which would explain the difference in saturated film thickness, 
concerns the balance of supply of matter through the creeping film and 
evaporation. Above the lambda point, where no film creep occurs, a 
thick film would evaporate faster than it could be formed from the bulk 
liquid or vapour. A question also arises as to the extent to which such a 
balance could determine the thickness below the lambda point. Although 
the contact of the mobile film is present, giving the observed thick film, 
it is conceivable according to these ideas that different conditions of heat 
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influx might lead to different values of the film thickness. If this were 
the case, the film would only be well defined in a dynamic sense in terms of 
its transport properties ; the static aspects such as thickness would have 
little significance. 

In evaluating the heat influx in our experiment, it should be remembered 
that the specimen tube is surrounded by liquid helium and liquid nitrogen 
in silvered dewars with narrow slits. Black paper was placed over the 
slits of the outer dewar to prevent radiation in the visible region from 
entering. This exclusion of visible radiation is to be contrasted with the 
experiments of Burge and Jackson and also Atkins. Accordingly, we 
have investigated the effect of shining light on the specimen. 


Haperiments with Radiation Influa 

For this experiment a 10w bulb was placed 20cm from the outer 
dewar in a galvanometer lamp housing. The lens was replaced by a 
shutter and a variable aperture. Filters could be placed in the light path. 

Our intention was to use a source of radiation which approximated the 

sources used in visual observations on the film flow (Atkins 1949) and 

optical investigations such as those of Burge and Jackson. Of course 
we have not reproduced exactly the optical system used by them and 
it is quite certain that the amount of radiation absorbed between the 
source and specimen is somewhat different. For this reason we shall 
not describe our results in quantitative detail, but rather concern our- 
selves with the general magnitude of the effect. 

In the first observations, a sheet of heat-absorbing glass and a heavy 
green filter were placed in the light path. With an aperture of 1-3 cm, the 
effect of the source was to evaporate the film to 40° of its original value ; 
the evaporation was complete in 1} minutes. It is important to demon- 
strate that the deflection of the balance which results from the light source 
is really due to an evaporation of the film and not tosome movement of the 
gas around the specimen resulting from local heating. This we have 
demonstrated in three ways. The radiation has been directed on to the 
counterweight ; the latter has a very small surface area so that evaporation 
of the film would be unobservable but motion of the ambient gas could 
deflect the balance appreciably. No deflection was observed. It may be 
argued that the different geometrical form of specimen and counterweight 
reduces the validity of this test. | Accordingly, we have made the following 
checks on the specimen. Having observed the evaporation of the 
saturated film at 1-80°K (vapour pressure=12-3 mm Hg), the temperature 
was altered to 2-02°K (vapour pressure= 24-8 mm Hg) and observations 
made on the unsaturated film at 50°% saturation. The pressure in the 
specimen chamber for 50% saturation at 2-02°K is the same as for complete 
saturation at 1-80°k, so that any deflection of the balance due to gas 
movement would be the same in both cases. However, the effect of 
film evaporation would be quite different because the saturated case 
involves 150 layers and the unsaturated 3 or 4. We indeed established 
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that the deflection due to radiation in the unsaturated case corresponded 
to less than 5 layers evaporated, hence the effect of gas motion is negligible. 
The only objection that can be made against this observation is that the 
presence of liquid might be of importance. This is not valid since we have 
observed that the same radiation produces an evaporation of about 3 layers 
in 13 above the lambda point in the presence of bulk liquid. 

The evaporation can be reduced by decreasing the aperture. 
Exchanging the green filter for a yellow filter resulted in little change, 
indicating that most of the energy is entering via the infra-red. The 
amount evaporated can be reduced by extra sheets of heat-absorbing 
glass. At an aperture of 1-3cm, one sheet of glass resulted in 60% 
evaporation, two sheets 50%, and three sheets 37%. 

In giving only the general magnitude of the effect, we hope to emphasize 
the fact that the actual values recorded by us have little significance except 
to the extent that the light source approximates that commonly used in 
liquid helium experiments.* A quantitative study of the effect would 
require some estimate of the energy impinging on the specimen itself and 
this we have not done. The evaporation from a specimen due to incident 
radiation would also depend on the nature of the contact between the 
specimen and the liquid. This contact may well be better in the case of 
Atkins’ experiment and also in Burge and Jackson’s than in ours, 
Accordingly, the effect of incident radiation would be less in their case. 

Burge and Jackson have considered the effect of radiation in their 
earliest experiment. One observation made by them involved turning off 
their light source for some time and then turning it on again to see if the 
film was any thicker. No increase in film thickness was detected. Two 
remarks are necessary on this observation. In the first place, we have 
observed the evaporation to be very rapid, the bulk of the film being 
evaporated in the first 30 seconds. More important than this, the 
observation (and others in the paper) obviously concerns only heat input 
by the light source and has no bearing on the more general question of 
whether differences in the heat influx in the experiments of various 
workers lead to disparities in values of film thickness. Such a possibility 
is reduced to a conjecture ; the sensitivity to heat influx is evident from 
our radiation experiments, but they do not prove the point. 

Burge and Jackson have drawn attention to the very rapid rise in the 
film thickness at the lambda point, particularly as a basis for discussing 
residual heat influx. They point out that if the thickness is governed by 
the balance of flow onto the specimen through the mobile film and 
reduction due to evaporation, then the slope of this thickness-temperature 
curve should change with heat influx. They observed no change as the 
result of turning off their light source. It is interesting that our observa- 
tion on the slope of this line agrees with that of Burge and Jackson. 
Such an observation appears to contradict the idea that: residual heat 


* The source used can have little effect on the transfer rate (Bowers and 
Mendelssohn 1950). 
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influx is responsible for the disparity between our results. We have found 
some difficulty in investigating the region close to the lambda point, but 
we hope to investigate this point further in future experiments, since 
Burge and Jackson’s argument is an important one. 

In our previous paper (Bowers 1953 b) we discussed the observed 
mobility of saturated and unsaturated films (Bowers, Brewer and Mendels- 
sohn 1951, Long and Meyer 1952 a, 1952 b) in terms of the film thickness 
involved. A further discussion of this question will not be undertaken here 
because it must first be established that there is some definite saturated 
film thickness which is independent of experimental conditions. 

In concluding this discussion, some remarks on the method used in this 
investigation are appropriate. This investigation has demonstrated that 
the gravimetric method may be applied to the problem of the saturated 
He II film. Much could be done to improve the sensitivity by providing 
space for a larger specimen. Furthermore, other surfaces may be | 
investigated. Unfortunately, the substrate of greatest interest for 
comparison with other work, i.e., glass, is the most difficult. The problem 
is that of shaping glass sheet to provide a large enough surface area with a 
small total weight. 
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ABSTRACT 

The possibility is considered that the 13-70 mev (2—) state of *°Ne 
is 7’—0 in which case the electric dipole transition to the 1-63 Mev (2-+-) 
state of 2°Ne is forbidden by the isotopic spin rule. This transition is 
sought and is shown to have a radiative width of less than 2 ev 
corresponding to a (2J-+1) | M |? value of less than 0-007. From this 
we may deduce that the contamination of this state with 7’=1 is less 
than 3% in intensity if we neglect the possible diminution of the transition — 
matrix element due to the high level density. 


§1. INTRODUCTION 

THIS paper continues the series (I—Wilkinson and Jones 1953, 
II—Wilkinson 1953 a, I1I—Clegg and Wilkinson 1953) in which an 
attempt is being made, by investigating the validity of the isotopic spin 
selection rules, to assess the purity of isotopic spin states in a variety of 
light nuclei with various degrees of excitation. As I, II and IIT the 
present paper concerns itself with the rule that electric dipole radiation 
is forbidden between states of the same isotopic spin in self-conjugate 
nuclei (Radicati 1952, Gell-Mann and Telegdi 1953). 


§2. States oF °Ne 

Many states are known in *°Ne above the excitation of 12-87 Mev 
produced by adding a proton to *F. Several of these states decay by 
emission of alpha-particles to states of 1®O at 6 or 7 Mev excitation ; 
we concern ourselves principally with that of 13-70 Mev formed by 
protons of 874 kev ; it is (2—) (Ajzenberg and Lauritsen 1952, Seed and 
French 1952), The total width of this level is 5-2 kev made up of a 
proton width of 1-1 kev and an alpha-particle width of 4:1 kev (Chao 
1950) or vice versa. It matters little which width we take as which ; 
we adopt the first of the alternatives since it is that usually chosen. 
This corresponds with a conventional interaction radius of 5-7 x 10-33 em, 
to reduced alpha-particle widths of about 0-8, 0-8 and 2% of the Wigner 
limit (Teichmann and Wigner 1952) for the three alpha-particle groups 
involved. It therefore seems probable that this state is of 7=0.+ 
We cannot be certain of this assignment since the density of known 
states in °F at the excitation of 3-4 Mev that corresponds to the region 


* Communicated by the Authors. 
{Compare reduced widths of 0-52, 0-93, 0-03, 0-12 and 0-016% of the 


Wigner limit for the indubitable 7'=0 states of 2°Ne at 6-7 to 7-9 Mev 
(Cameron 1953). 
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of Ne of present interest (the first 71 state of 2°Ne is expected at 

about 10 Mev) is high (about 6 per Mev) and comparable with that observed 
for ?°Ne (about 12 per mev). The possibility that we have guessed 
wrongly will be considered later. We may also consider the state in 2°Ne 
at 13-44 Mev formed by protons of 598 kev; according to unpublished 
work by J. Seed and G. Dearnaley this state is also probably (2—) and, by 
the same argument as that presented above, may well also be 7'=0. 
(The widths for this state are 0-04 and 37 kev.) 

The first excited state in ?°Ne is at 1-63 Mev and is certainly of even 
parity from the ft value of the B-decay of »°F (Alburger 1952 and 
Bromley, Bruner and Fulbright 1953) and is expected to be (2+) 
(Scharff-Goldhaber 1953); this expectation has been confirmed by 
Seed (1953). There is therefore the possibility of electric dipole radiation 
of about 12 Mev between the (2—) and (2+) states. But if we are 
correct in our assignment of 7’'=0 to the (2—) states these transitions 
should be isotopic spin forbidden. They have indeed long been known 
to be feeble compared with the well-known capture radiation from the 
13-51 Mev state formed by protons of 669 key (Devons and Hereward 
1948) ; J. KE. Sanders (1952) has reported that the ratio of their intensity 
to that of alpha-particles leading to excited states of 1°O is less than 5% 
of the same ratio for the 13-51 Mev state. 


§ 3. EXPERIMENTAL INVESTIGATION 


We have carried out a survey of the capture radiation from several 
states of 2°Ne including those at present under consideration. We have 
used Nal(Tl) crystals for this work and have sought the expected 
12 Mev radiation directly against the background of the much greater 
flux of 6 and 7™Mev radiation from the states of 16O reached by 
alpha-particle emission. The crystal was at 90° to the proton beam ; 
the expected angular distribution for the 12 mev radiation is 1+;4 cos? 6. 
This work will be reported in detail elsewhere but we may say here that 
the well-known uncertainties associated with the interpretation of a 
high-energy ‘tail’ in the crystal pulse spectrum were reduced by 
irradiating the crystal under identical conditions with the 7:48 Mev 
gamma-rays from the reaction °Be (py)’°B and those of 8-06 Mev given 
by C(py)“N. The 1°B and “N distributions, which were similar in 
form, were used to construct the expected tail due to the 6-1 and 7:1 
radiations from 1F (p«)!®O; which tail, correctly normalized, was 
compared with the distribution due to the bombardment of fluorine by 
protons when any excess at high energy in the latter distribution would 
represent the contribution from gamma-rays of energy greater than 
7mev. There is no significant excess and upper limits of 0-05 and 
0:04°% maybe placed on the intensity of the expected 12 Mev radiation 
relative to that of the 6 and 7mev gamma-rays for the 13-70 and 
13-44 ev states respectively (using the data of Devons and Hereward 
(1948), Rae, Rutherglen and Smith (1950), Carver and Wilkinson (1951), 
and our own work for the relative abundances from the 13-51 Mey state). 
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$4, INTERPRETATION 


If we continue our assumption that [,>J), our result implies that 
Lee and 15 ev for the electric dipole radiation of 12 mev from the 
13:70 and 13:44 Mev states respectively. This corresponds to values 
of (2J-+1) | M |?<0-007 and 0-05 for these two transitions in terms of 
the estimate of the single particle model (Weisskopf 1951). Now for 
allowed electric dipole transitions in light elements (2J+-1) | M |?-~0-2 
with a rather small spread (Wilkinson 1953 b) and so the first of the 
present transitions appears to be discouraged by a factor of about 30 or 
more, while little may be said about possible discouragement of the second, — 

It may be argued, however, that when the level spacing becomes small 
radiative widths should become smaller than those given by the 
single-particle model (see, for example, Blatt and Weisskopf (1952)) ; 
this effect has been noticed by Kinsey (1953) in a survey of experimental 
material throughout the periodic table. Kinsey suggests | 1 |?~0-1 D 
where D is the level spacing (in Mev) for states of the same spin and parity. 
Now at the present excitation in ?°Ne D~300 kev so Kinsey’s suggestion 
would yield for the (2—) state in question (2J+1) | |?~0-15 in close 
accord with the earlier suggestion.* It is therefore doubtful whether 
‘the effect of level density is yet of great importance; (2J+1) | M / 
values similar to those for levels of a few Mev excitation are found for | 
transitions from states at 16-1 Mev in #C, 13-1 Mev in 1%O, 13-51 Mev 
in 2Ne, 12:6, 13-0 and 13-0mev in *°Si where the level spacings 
(irrespective of spin and parity) are respectively about 500, 200, 100, 30, 
20 and 20 kev. Conversely no cases are established with certainty in 
light nuclei of high level density of dipole radiative widths that fall 
greatly below the single-particle figure (though two possibilities of 
magnetic dipole transitions with (2J-+-1) | M |?<0-1 occur in *Ne—from 
the 13-19 Mev state). 

If, then, we may interpret the above ‘ discouragement factor’ as due 
to the isotopic spin rule we would say that the contamination of the 
predominantly 7=0 state at 13-70 Mev with 7’'=1 is less than 0-03 in 
intensity; in the notation of Radicati (1953) «92(1)<0-03.f We must 
stress, however, that this conclusion is rather more tentative than those 
based on the behaviour of more widely-spaced levels at lower excitation 
owing to our uncertainty of the exact importance of the increase in level 
density for the uninhibited transition strengths. It is perhaps a little 
- surprising that the contamination should be as low as 0-03 since, as we 


* It must be emphasized that there is little or no justification for the factor 
(2J-+-1) used in these comparisons ; its use is purely one of convenience in 
reducing and tabulating the experimental data, since it is usually (2J-+1)I, 
rather than I, that emerges from the experiments. 

j If the choice ,>T, turns out to be the right one for this level, the 
values of Ty, | W|* and a, (1) must be diminished by a factor of four, 
This makes little difference to the argument. If, however, [,)>T,, for the 
13:44 Mev state we find, for the forbidden electric dipole tzansition, 
Ty<0-02 ev, (2J-+1) | M |?<0-0001, a ?(1)<0-0005 ; these values appear so 

. low as to be inadmissible. 
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have noted above, D~300 kev, which is of the same order as the probable 
matrix element of the Coulomb perturbation that is responsible for 
mixing states of different isotopic spin. We may compare the situation 
for the 13-09 mev state in 10 where, with a value of D of about two or 
three times greater than in Ne, «,2(0)>0-03 (Wilkinson 1953 c). 

We may, as noted in §1, have guessed wrongly that 7=0 for the 
13-70 Mev state. If in fact 7=1 the electric dipole transition would 
be allowed by the isotopic spin rule and we must regard the low value 
of (2J+-1)| M |? as due either to the effect of level density or to an 
unusually low value for the matrix element or to both, and this may 
appear unlikely, although it is certainly possible. For 7’=1 the emission 
of alpha-particles would be forbidden, and we would say, following the 
data presented in §1, «,?(0)>0-04* with the strong possibility that 
a more realistic estimate would be considerably larger since this 
figure is based on the Wigner limit and it is rather unlikely that 
a contaminating level at so high an excitation in ?°Ne should possess 
single-particle alpha-particle characteristics.+ 
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SUMMARY 


A comparison is made between the nuclear disintegrations observed 
when 140 Mev protons pass through photographic emulsions and the 
nature of the disintegrations expected on the Goldberger model (1948). 
Satisfactory agreement is found for all the features of the disintegrations 
examined. 


§ 1. INTRODUCTION 
In a previous paper from this laboratory (Lees et al. 1953), referred to as I, 
the cross sections for the production of nuclear disintegrations by 130 Mev 
protons in the heavy (silver and bromine) and light (carbon, nitrogen and 
oxygen) constituents of photographic emulsions were obtained. The 
results were shown to be in agreement with the ‘optical model’ of Fernbach, 
Serber and Taylor (1949). 

The* optical model ’ assumes that at high energies a nuclear interaction 
can be produced by the collision of the incident nucleon with a single 
nucleon inside the nucleus, and merely predicts the probability for the 
occurrence of this collision. After the first collision, however, both 
nucleons may have sufficient energy to undergo further collisions inside 
the nucleus (i.e. initiate a nucleon cascade), and the knocked-on nucleons 
can reach the surface of the nucleus and escape, or lose sufficient energy to 
be absorbed in the nucleus leaving it in an excited state. Charged and 
uncharged particles are subsequently evaporated from the excited 
nucleus (Weisskopf 1937). 

A detailed theoretical description of the cascade process has not so far 
proved feasible. Goldberger (1948) has, however, suggested that numeri- 
cal results may be predicted by considering the passage of a sufficient 
number of nucleons through the nucleus and allowing the position and 
nature of the collisions to be determined by the laws of chance. This will 
be referred to as the ‘ Goldberger Method ’. 

Bernardini et al. (1952) have shown that the Goldberger model of 
successive single nucleon-nucleon collisions inside the nucleus is satis- 
factory for describing the stars produced by 400 Mev protons in photo- 
graphic emulsions. It may not be as reliable for primary particles of 


* Communicated by Professor P. I. Dee, F.R.S. 
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lower energy, since the cross sections for free nucleon-nucleon collisions 
increase as the energy of the primary particles becomes smaller, and at the 
same time their de Broglie wavelengths increase. It is, therefore, of 
interest to determine the range of validity of the single nucleon-nucleon 
picture. In this paper an account is given of a calculation carried out by 
the Goldberger method for the stars produced by 140 Mev protons in 
nuclear emulsions. The calculated data will be compared with the 
experimental results obtained using G5 emulsions which were exposed to 
146 Mev protons accelerated by the Harwell cyclotron. 


§ 2. METHOD oF CALCULATION 


The calculations were carried out for a heavy nucleus (Z=41, A=94), 
and a light nucleus (Z=7, 4=14). The nuclear radius was taken as 
1-42 10°%A'/3 cm. One hundred primary protons were followed in 
each case. 

A Fermi gas model was used for the nucleons inside the nucleus ; their 
maximum energy was taken as 22Mmev. The binding energy of the 
nucleons was taken as 8 Mev, making the depth of the potential well 
equal to 30mMev. It was assumed that the incident protons moving 
with a kinetic energy of 140 Mev gained 30 Mev on entering the nucleus. 

In making the calculation it was assumed that the primary proton 
could initiate a nucleon cascade inside the nucleus, and that the cascade 
would develop through a series of single nucleon-nucleon collisions. 
Allowance was made for the following features in each encounter :— 

(a) collision with either a neutron or proton, 

(6) variation of the free n—p and p-—p scattering cross sections with 
energy (Christian 1952, Pickavance and Cassels 1952), 

(c) variation of the energy of struck particle and angle between the 
colliding particles inside the nucleus, 

(d) variation of the angular distribution of n—p and p-—p scattering in 
the centre of momentum system with energy (Christian loc. cit., Randle 
et al. 1952). 

The detailed procedure for carrying out the calculation will not be 
described, as it was essentially the same as that used by Bernardini et al. 
(loc. cit.). 

The individual nucleons set in motion during the cascade were followed 
until they reached the nuclear surface, or were absorbed. Reflection at 
the nuclear surface was neglected. A value of 6 Mev, which is smaller 
than the accepted value of 8-9 Mev, was used for the Coulomb potential 
barrier for protons in the heavy nucleus. This figure was adopted so as 
to allow for some Gamow penetration of the barrier. Therefore, if after 
a collision the energy of a proton was less than 36 mev (inside the heavy 
nucleus), or 30 Mev in the case of a neutron, it was assumed that the 
particle was absorbed leaving the nucleus in an excited state. Thirty 
Mev were subtracted from the energy of each particle as it left the nucleus. 
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The excitation energy of the nucleus following the cascade is given by 
N 
U=[140— 2 #,—8(N—1)] Mev, . = == 
i=1 


where E, is the energy of a knock-on nucleon (after escaping from the 
nucleus), and N is the total number of knock-on nucleons emerging from 
the nucleus. 

The calculations were not extended to nucleons with energies below 
60 Mev inside the light nucleus (see § 5). 
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Mean free paths for nucleons in nuclear matter ; fo) points derived from cascade 
calculation ; value expected from free nucleon-nucleon scattering 
without application of Pauli exclusion principle ; — — — — value modified 
by Goldberger’s formula to allow for the effect of the exclusion principle. 


It was assumed that only single nucleon-nucleon collisions took place 
at all energies, and, therefore, at all energies the free particle cross sections 
were used. If in any encounter either particle was calculated to possess 
an energy less than 22 Mev after the collision, then it was assumed that 
the collision was forbidden by the Pauli Exclusion Principle. It is of 
interest to note the effect of the exclusion principle as a function of 
energy. In fig. 1 the mean free path for nucleons in nuclear matter 
resulting from the calculations is shown as a function of energy. For 
comparison the uncorrected mean free path of nucleons in nuclear matter 
and the same path length corrected using an analytical expression given 
by Goldberger (loc. cit. eqn. (14)) are shown. It can be seen that the 
mean free paths obtained using this expression are rather lower than those 
obtained by us. This difference is to be expected since Goldberger 
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ieee isotropic scattering in the centre of momentum system when 

eriving we formula, whereas the rise observed experimentally near 0° 
and 180 m n-p scattering tends to increase the effect of the Pauli 
Exclusion Principle. 


§ 3. THE Cross Sections ror Nuciear INTERACTIONS 


In table 1 the cross sections for nuclear interactions in light and heavy 
nuclei derived by the Goldberger method are compared with those 
obtained in Paper I. It can be seen that the overall agreement is 
satisfactory. 


Table 1 

Heavy Nucleus Light Nucleus 
No. of primary 
protons 100 100 
No. of nuclear 
interactions 81 59 Goldberger 
Cross section for Hu De. 
nuclear inter- 
action 4 (104 +15) x 10-26 em? | (2143) x 10-26 cm? 
Cross section for 
nuclear inter- 
action—Experi- 
mental (1282-16) ><310=2° cm? | (162-7) x 10-** em? 
Cross section for Paper I 


nuclear inter- 
action—Perry’s 
Method (1952) (110 +5) 10-26 cm? | (24+2) x 10-26 cm? 


§ 4, COMPARISON OF THE CALCULATED AND EXPERIMENTAL DATA FOR 
PRotToNS EMITTED WITH ENERGIES GREATER THAN 30 MeV 


The model of successive single nucleon—nucleon collisions inside the 
nucleus is most likely to hold for high energy particles, whose mean free 
paths for such collisions are considerably greater than the average 
inter-nucleon distance. This is considered to be the case for nucleons 
of energy greater than 60 Mey inside the nucleus (Peaslee 1952). Further- 
more some of the charged particles leaving the nucleus with energies | 
less than 30 Mev are emitted during the evaporation process. The number 
of evaporation protons with energies above 30 Mev is negligible. Therefore 
in this section we will confine our comparison to protons emitted with 
energies greater than 30 Mev, and will extend our analysis to energies 
below 30 Mev in § 5. 

The results obtained by the Goldberger method for heavy and light 
nuclei were combined in the proportion of the geometrical area presented 
by each to protons in Ilford G5 emulsions, that is approximately 3 to 1. 
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The experimental results were based on stars obtained by area scanning 
and were confined to those stars whose primary protons had energies 
greater than 130 Mev, their mean energy being close to 140 Mev. The 
energies of the primary and secondary protons were determined by grain 
counting. Some of the plates were underdeveloped so that energies as 
low as 30 Mev could be determined by this technique. 


4.1. Energy Distribution of the High Energy Protons 


In table 2 a comparison is made between the calculated and experi- 
mental distributions of protons emitted with energies greater than 
30 Mev. The data has been normalized in each case to represent the 
number of protons emitted per star. In Paper I it was pointed out 
that Op stars are difficult to detect in ‘ area ’ scanning, and that primary 
protons which lose less than 25 Mev in a nuclear encounter and are 
scattered through small angles are difficult to distinguish from elastically 
scattered protons. The effective number of stars in the calculated data 
was therefore made to correspond to experimental conditions by rejecting 
from the calculated results all the Op stars, and the 1p stars in the above 
category. 

Table 2 


Energy in Mev 30-50 50-70 70-90 90-110 >110 


No. of protons 
per star—Expe- 
rimental 0-23 +0-04 | 0-13 +0-02 | 0-07 +0-02 | 0-11 +0-03 | 0-05 +0-02 


No. of protons 
per star—Calcu- 
lated 0-22 +0-04 | 0-16 +0-04 | 0-21 +0-04 | 0-12 +0-03 | 0-08 +0-03 


4.2. Angular Distribution of the High Energy Protons 

In making the cascade calculations a two-dimensional geometry was 
used. Thus there is no exact correspondence between the calculated 
and the observed angular distributions. For a sufficiently large number 
of stars, however, the calculated angular distribution should approximate 
to the distribution of the actual angle in space between the primary 
proton and the secondary protons. The two distributions are compared 
in table 3, where the intervals given represent angles between the direc- 
tions of the primary and secondary protons. The calculated results are 


normalized to the total measured number of protons with energies 
greater than 30 Mev. 


4.3. Discussion 


It can be seen from tables 2 and 3 that the measured energy and 
angular distributions of the high energy protons are in satisfactory 
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agreement with the distributions calculated on the assumption that the 
primary proton initiates a nucleon cascade. The only point at which a 
large discrepancy appears between the experimental and calculated 
results is in the 70-90 Mev energy interval in-table 2. This difference 
could possibly be accounted for by statistical fluctuations, since both 
values deviate from the general trend of figures as a function of energy, 
and the deviations are in opposite senses. 


Table 3 


Angle in degrees | 0-20 | 20-40 | 40-60 | 60-80 | 80-100 | 100-120 | +120 


No. of protons— 
Experimental 


| No. of protons— 
Calculated 


§ 5. COMPARISON OF THE CALCULATED AND EXPERIMENTAL Data FoR 
THE CHARGED PARTICLES EMITTED WITH ENERGIES LESS THAN 
30 Mev FROM STARS PRODUCED IN SILVER AND BROMINE 


For nucleons of low energies (60 Mev inside the nucleus) the compound 
nucleus picture is believed to hold (Peaslee loc. cit.). In this model the 
incident nucleon is assumed to share its energy quickly with the rest of 
the nucleons in the nucleus giving rise to a compound nucleus which 
subsequently decays. If this model were entirely correct one would 
expect very few knock-on protons of energies less than 30 Mev to emerge 
from the nucleus. The model of single nucleon-nucleon collisions is also 
more likely to break down at low energies where simultaneous collisions 
of more than two nucleons may take place. In making our calculations 
by the Goldberger method, however, the nucleons were followed through 
the nucleus assuming single nucleon-nucleon collisions, until the particles 
either left the nucleus or were absorbed. A large proportion of the 
protons with energies less than 30 Mev then come out of the nucleus as 
their apparent mean free path inside was increased due to the action 
of the Pauli Exclusion Principle (see fig. 1). It is therefore interesting 
to compare the results of the calculation made on the assumption of a 
nucleon cascade with the experimental data for protons emitted with 
energies below 30 Mev. 

In this section the results of the calculations on the heavy nucleus 
will be compared with the experimental data. A fairly reliable evapora- 
tion theory is available for silver and bromine (Le Couteur 1950, 1952) ; 
such a theory would not be appropriate for light nuclei because of the 
limited number of nucleons, and the ease with which these nuclei disin- 
tegrate completely into «-particles and protons at comparatively low 
excitation energies. 
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In Paper I it was shown that 95%, of the 1, 2, 3p stars and 70% of the 
4p stars were produced in silver and bromine. Therefore, it is reasonable 
to compare the results of the calculations for the heavy nucleus with the 
sum of the experimental results for 1, 2, 3 and 4p stars. Protons from 
stars in this category having one or more «-particle of energy less than 
9 Mev were not included in the experimental data in §§ 5.2 and 5.3 as 
these stars were probably produced in the light nuclei (Paper I). 

All protons of energy less than 30 Mev will be called black. 


5.1. The Size-Frequency Distribution of Stars in Silver and Bromine 


For each interaction the Goldberger method gives the number of 
knock-on protons and the excitation energy of the residual nucleus. 
The calculated excitation energies varied from 2 to 148 Mev, the mean 
value being 44 Mev. Le Couteur (1950, 1952) has given a relation between | 


Table 4 


Calculated Calculated cross Experimental cross 
number of stars | section «1076 cm? | section x 107° cm? 


Type of Star 


14 
29-5 
39 
17 
4-5 


Totals 


the mean number of charged evaporation prongs and excitation energy. 
From this relationship the number of charged particles emitted during 
evaporation can be calculated for each star separately. The total number 
of charged prongs per interaction can thus be obtained. The resulting 
size-frequency distribution is given in table 4, where a comparison is 
made with the experimental results derived in Paper I. 


5.2. The Energy Distribution of Protons emitted in Silver and Bromine 

The experimental energy distribution of protons from 1, 2, 3 and 4p 
stars is shown in fig. 2; this was determined from the ranges of the 
protons ending in Ilford G5 emulsions. It is based on measurements 
on 388 stars ; the total number of black protons emitted in these stars 
after applying a suitable geometrical correction for those leaving the 
emulsion was 420. For comparison the calculated data for the heavy 
nucleus was normalized to 388 stars ; a correction was made for the Op 
and Ip stars to allow for those that would not be found by area scanning 
as discussed in § 4.1. 
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On the basis of our calculations by the Goldberger method we would 
expect 180 protons to be emitted from 388 stars in evaporation processes. 
In obtaining this result, allowance was made for the variation of the 
excitation energy, U, as discussed in § 5.1. The number of evaporation 
protons expected using the above model is, therefore, not sufficient to 
account for the 420 tracks observed experimentally. However, the 
Goldberger model also predicts that an additional 190 protons of energy 
below 30 Mev would be emitted from the heavy nuclei in direct collision 
processes. The total calculated number emitted is, therefore, in satis- 
factory agreement with the experimental number. 


Fig. 2 
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In fig. 2 the shaded area represents the calculated number of knock-on 
protons emitted in the energy range 6 to 30 Mev. The agreement with 
the experimental energy distribution is reasonable for energies greater 
than 20mev. Below 20mev the effect of the knock-on protons is 
completely masked by the evaporation protons. ; . 

According to simple evaporation theory (Weisskopf /oc.cit.) the 
probability for the emission of a particle of energy between # and 
E-+dE is Wiese 


E—V 
P(E )\dE= 73 exp [GP law tn ad 
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where V is the height of the Coulomb barrier and 7’ is the nuclear tem- 
perature of the residual nucleus. Assuming a Fermi gas model we have 


UapTt - P 1 ee ee 


where, according to Le Couteur, 8 is ~7-6 Mev~? for a mean silver and 
bromine nucleus. In order to compare the calculated with the experi- 
mental results values must be assumed for V and 7’. In eqn. (2) 7 
represents the temperature of the nucleus after emitting a proton. The 
excitation energy varied between 2 and 148 mev, the mean value being 
44mev. At low excitation energies neutron emission predominates so 
that most of the protons observed are emitted at an average excitation 
energy greater than 44 Mev. However, when the total number of evapor- 
ation protons was calculated, allowance was made for the variation of 
U as mentioned above. Therefore, the uncertainty in U affects the form 
of the evaporation spectrum only, but not the total number of evaporation 
particles emitted. Furthermore, 7 is proportional to »/U so that the 
uncertainty in the excitation energy of the residual nucleus does not 
seriously affect the shape of the energy spectrum of the protons resulting 
from evaporation. A value of 2-4 Mev was taken for 7’ corresponding 
to an excitation energy of the residual nucleus of 44 Mev. 

The experimental data presented in fig. 2 suggests that the effective 
Coulomb barrier is about 4 Mev. The potential barrier for a mean silver 
and bromine nucleus (unexcited) is 8-9 Mev for protons. However, the 
effects of barrier penetrability and thermal expansion of the excited 
nucleus could reduce the effective height of the barrier to ~4 Mev. The 
resulting curve for 180 evaporation protons with 7’=2:4 and V=4 Mev 
is shown in fig. 2; it is also shown added to the calculated knock-on 
spectrum. The agreement of the combined spectrum with the experi- 
mental spectrum is satisfactory, except that some protons of energies 
less than 4 Mev are found experimentally. Some of these may be from 
1, 2, 3 and 4p stars produced in light nuclei, from which no «-particles 
of energy less than 9 Mev are emitted. Some of the particles may also 
be ‘decay protons’ from heavy nuclei (Le Couteur 1950, Fujimoto and 
Yamaguchi 1949). 

The results obtained by the Goldberger method are, therefore, adequate 
to account for the numbers and energy distribution of protons below 
30 Mev emitted from stars produced by 140 Mev protons. 

It should be pointed out, however, that processes different from those 
assumed in our calculations could also account for the experimental 
energy spectrum below 20 Mev. If it were assumed that all nucleons 
with energies below 50 Mev (inside the nucleus) were absorbed, then the 
average nuclear excitation energy would be raised from 44 to 66 Mev. 
On using Le Couteur’s theory, as described in § 5.1, the total number of 
protons expected to be evaporated with energies between 0 and 20 mev 
would be 310 compared with 3840 protons observed experimentally. 
The energy spectrum calculated on this assumption is shown in figs dc: 
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If it were assumed that all nucleons of energy less than 60 Mev inside 
the nucleus were absorbed, the average excitation energy would be 
increased still further but the form of the calculated evaporation spectrum 
would not be substantially different from that shown in fig. 3, and would 
not be in agreement with the number of protons in the energy range 
20-30 Mev observed experimentally.* It therefore seems probable that 
some knock-on protons of energies between 20 and 30 Mev escape from 
the nucleus. Additional evidence for this view will be presented in the 
next section. 

Fig. 3 
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5.3. Angular Distribution of Protons of Energies less than 30 Mev emitted 
from Stars in Silver and Bromine 


The protons emitted during the evaporation of a heavy nucleus should 
be ejected isotropically relative to the recoiling nucleus. The velocity 
of the evaporation protons is ~0-1 ¢ whereas in our case the mean velocity 
of the recoiling nucleus is ~0-005c. Therefore the evaporation protons 
should have an approximately isotropic distribution in the laboratory 
system. The experimental angular distribution for protons of energy 
less than 30 Mev from 1, 2, 3, and 4p stars shows a marked anisotropic 
distribution, the ratio of the number of black protons in the forward 


* Data from 900 stars was used to determine the mean number of protons 
in the energy range 20-30 Mev. The energies were determined from grain 
counts in the underdeveloped plates and from the ranges of the protons ending in 
theemulsion. The mean ordinate normalized to 388 stars in figs. 2 and 3 was 17 +4. 
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direction to the number in the backward direction being 1-9+-0-2. The 
calculated forward to backward ratio obtained by combining the calcu- 
lated angular distribution for black knock-on protons with an isotropic 
angular distribution for the evaporation protons is 1-9+-0-5. This value 
is in good agreement with the experimental figure and confirms that some 
knock-on protons of energy less than 30 Mev escape from the nucleus. 

It was shown in the previous section that an evaporation process 
could not easily account for the production of protons with energies in 
the region from 20 to 30 Mev. The forward to backward ratio for the 
protons in this energy range offers confirmation for this conclusion ; 
20 protons were observed travelling in the forward direction, and 3 in 
the backward. 

It is of interest to extend these observations to protons of energies 
less than 20 Mev. The Goldberger model gives a forward to backward 
ratio of 1-7-+0-5 for protons of energies less than 20 Mev. *The experi- 
mental forward to backward ratios for protons ending in the emulsion 
are 1:5+0-2 and 1:4+0-2 for protons in the energy ranges 0-20 Mev 
and 0-10 Mev respectively. However, some of the protons included in 
these figures could have come from stars produced in the light nuclei from 
which no «-particle of energy less than 9 Mev was emitted. Protons from 
identifiable stars produced in the light nuclei show a marked anisotropic 
distribution, so that it is not possible to conclude with certainty that 
knock-on protons of energies less than 20 Mev are emitted in the stars 
produced in silver and bromine. 


5.4. The Ratio of the Number of «-particles to the Total Number of Particles 
emitted in the Stars produced in Silver and Bromine 
The experimental ratios for the number of «-particles to the total 
number of prongs in a star as a function of star size were presented in 
Paper I (fig. 2), and the values for 1, 2 and 3p stars are reproduced in 
table 5. This ratio will be referred to as the «/total ratio. 


Table 5 


a/Total 


Type of star 


Calculated for Ag and Br Experimental 


Ip 0-04 40-02 0-05 40-02 
2p 0-08 0-02 0-07 40-01 
3p 0-10 +0-03 0-16-+9-02 


According to Le Couteur (1950), the ratio of the number of a-particles 
to the total number of charged evaporation particles emitted from silver 
and bromine should be 0-25. This figure was used in calculating the 
a/total ratios ; the total number of charged particles in each star was 
obtained as discussed in § 5.1. The results are presented in table 5. 
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The agreement for 1p and 2p stars is very good, and the low a/total 
ratio found experimentally again shows that an appreciable number of 
knock-on protons are emitted in the stars. If it were assumed, as in 
§ 5.2, that all nucleons of energy less than 50 Mev inside the nucleus 
might be absorbed, the increase in the excitation energy would approxi- 
mately double the calculated number of charged evaporation prongs 
emitted, thereby doubling the calculated «/total ratios. It can be seen 
that the experimental «/total ratios for 1p and 2p stars are in better 
agreement with the theory that some knock-on nucleons of energies less 
than 50 Mev (inside the nucleus) escape from the nucleus. The calculated 
a/total ratios for 3p stars given in table 5 are significantly lower than the 
experimental value for 3p stars. This is probably due to the fact that 
some of the 3p stars found experimentally were produced in light nuclei 
for which the measured ratio is ~0-5, rather than to the absorption of 
nucleons of energies less than 50 Mev inside the nucleus. 


5.5. Discussion 


From the number of protons occurring with energies greater than 
20 Mev as shown in fig. 2, and the experimental forward to backward 
ratio for protons in the energy range 20-30 Mev, discussed in § 5.3, it 
is concluded that some nucleons with these energies escape during the 
nucleon cascade. There is no definite evidence for the existence of 
knock-on protons with energies below 20 Mev, but the forward to back- 
ward ratio for protons less than 20 Mev produced in 1, 2, 3 and 4p stars, 
and the «/total ratios found experimentally for Ip and 2p stars suggest 
that some nucleons of energies less than 20 Mev may also escape during 
the nucleon cascade. 

This result suggests that nucleons of moderate energy may travel 
appreciable distances inside a nucleus before being absorbed. Evidence 
supporting this view has been recently presented by Feshbach, Porter 
and Weisskopf (1953). In making an analysis of the experimental data 
of Barschall and his co-workers (1952) on the total cross sections of heavy 
nuclei for the scattering of neutrons with energies between 0-1 and 
3 mev, Feshbach et al. have concluded that the incident neutrons penetrate 
a distance of about 2x10-!2cm into nuclear matter before they are 
incorporated into the compound nucleus. This distance may be compared 
with the figure of ~10-!cm for the mean free path of low energy 
nucleons shown in fig. 1. 

There is a possibility that sometimes a nucleon collides simultaneously 
with more than one nucleon inside the nucleus, leading to a sharing of 
energy between three or more particles. The simplifications and poor 
statistics inherent in the calculations do not allow us to exclude this 
possibility for nucleons of energies less than 60 Mev inside the nucleus, 
Though probably not entirely correct in detail, the model of a mucleonic 
cascade of single nucleon-nucleon collisions together with Le Couteur’s 
theory of the evaporation process appears to be adequate to account for 
the main properties of stars produced in silver and bromine. 
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CXLI. On the Annihilation of Positrons in Solids 
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ABSTRACT 


Using the method of delayed coincidences the lifetimes of positrons 
from ??Na have been determined in various materials, the 1:28 mev 
y-ray of Na announcing the emission of a positron. 

A lifetime of (1-6+-0-6 x 10-1° sec) has been deduced in aluminium, and 
comparative measurements show that it does not differ from this in other 
metals studied by us and by other writers. A second lifetime of 
(2-5-L0-3 x 10~® sec) involving (0-27-0-03) of the positrons has been found 
in distrene at room temperature. This lfetime has been found to increase 
by about 10% for 115°c rise in temperature. Three quantum emission 
from distrene has been found to be small. Experiments have also been 
carried out in plastic sulphur and magnetized iron. 


§ 1. INTRODUCTION 


Tue theory of positron annihilation in its early form (Heitler 1944) 
predicted a positron lifetime in lead of about 10~1° sec and lifetimes in 
other materials proportional to their electron density, assuming all Z 
electrons to be available. The experiments of De Benedetti and Richings 
(1951) showed that there might be little variation from metal to metal. It 
seemed therefore desirable to investigate both comparative and absolute 
positron lifetimes in a variety of materials. This would indicate how 
closely the theory is followed. 

The effect of a magnetic field in iron on positron lifetimes was studied. 
This work arose from consideration of the effects of spin on the annihilation 
process. ; 

The experiments of Deutsch (1951) on annihilation in gases in which 
positronium formation was established, suggested that complex time 
distributions might be found in the annihilation radiation from solids. 
Further it indicated the possibility that three-quantum annihilation might 
be more prevalent than the 1/370 of two-quantum processes calculated on 
the basis of free positron electron annihilation (cf. Ore and Powell 1949). 

Bell and Graham (1952) had reported a complex time distribution in 
certain plastics and a difference in behaviour between crystalline and fused 
quartz. The latter difference pointed to considerable influence of the gross 
On a pe cae Ee nds a li ea a ate SE a rae ae 
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structure on annihilation. From these considerations it was decided to 
investigate the lifetimes of positrons in distrene and their variation with 
its expansion, melting and previous treatment. A search for three 
quantum annihilation was also carried out. 

Lifetimes in the allotropic modifications of sulphur seemed to be of 
interest in the light of the reported behaviour of quartz. 

In the meanwhile De Benedetti and Richings (1952) and Bell and 
Graham (1953) have published further work in this field. 


§ 2, APPARATUS 


The lifetimes were measured by the method of delayed coincidences 
using a coincidence unit of resolving time t>=2 x 10~® sec. 

The principles of this unit are essentially the same as that of Bell, 
Graham and Petch (1952). A block diagram of the arrangement is shown 
in fig. 1 (ef. Lewis and Ferguson 1953). The scintillators were solutions 


Fig. 1 


| Poteet Mans eee Sees: 


Block diagram of apparatus. Full line shows basic apparatus. Dashed line 
shows later additions. C.1, 2—counters; A.1, 2, 3, 4—amplifiers; 
D.1, 2, 3—discriminators ; K—kicksorter ; 3-C—triple coincidence unit ; 
F—cathode follower ; G—gate ; d—1-25 psec delay line. 


of terphenyl in toluene contained in plane ended glass vessels 1 inch 
diameter by 1 inch long, surrounded by magnesium oxide. The photo- 
multipliers were 14-stage E.M.I. type 6262 selected for large gain and 
high photosensitivity. They were run at voltages between 2000 and 2500. 
To obtain a linear response it was necessary to take the side channel 
pulse from the llth dynode. Delay was introduced between the coin- 
cidence unit and counter by 952 Uniradio 31 coaxial cable which gives a 
delay of 5 x 10~® sec per metre. 


Annihilation of Positrons in Solids 1341 


The resolution curve, obtained by measuring the number of coincidences 
against delay introduced, shows a resolving time dependent particularly 
on the setting of the middle channel discriminator and times down to 
1x 10~° sec are obtainable, though the efficiency at this resolving time is 
very poor. The actual position of the resolution curve is sensitive to the 
settings of the side channel discriminators D.1 and D.2 and to the 
photomultiplier voltage. 


§ 3. Earty CoMPARISON EXPERIMENTS 


Using the apparatus described, at an earlier stage of development, 
comparison measurements were made of positron lifetimes in various 
materials. The ?*Na source used was enclosed in 0-005 inch aluminium. 
It was placed between samples of the material under test and lead, the 
standard throughout these measurements. The geometrical arrangement 
of source and counters was the same as that of De Benedetti and Richings 
(1952). 

The discriminator D.1 (see fig. 1) was set to exclude pulses from 
511 kev y-rays, and D.2 so as to minimize coincidences caused by scattered 
y-rays. This ensured that the 1-28 Mev y-ray of ?2Na was detected in 
counter 1, and an annihilation y-ray in counter 2. 

A resolution curve was obtained first with the sample then with lead, 
next the source. A typical set of data is shown graphically in fig. 2, 
comparing iron and lead. 

The results were obtained by the analysis of Bay (1951). No difference 
was found between the positron lifetimes in lead and aluminium, 
magnesium and iron. The error is estimated to be +5%x 10-11 sec. 
A centroid shift of about 2x 10-19 sec was found for wax, mica and water, 
indicating a lifetime of not less than this in these materials. With the 
resolving time used, measurements on the non-metals were not continued 
to sufficiently long delays to detect the complex behaviour reported for 
these non-metals by Bell and Graham (1953). 


§ 4. INVESTIGATION OF Lone PERIODS IN SOLIDS 


Tn the search for a long period in solids, first, the tails of the resolution 
curves for metals were examined at long delays but no evidence of the 
existence of any second period was found. 

The equality of the lifetimes in iron and the other metals indicates that 
the aligned spins of the d-shell electrons of the domain structure of the 
iron has little effect. To examine this matter further, the specimen, in the 
form of a hollow cylinder was magnetized. With a delay of 3x 10-° sec 
(see fig. 2) the number of coincidences was measured alternately with and 
without an induction of some 7000 gauss per cm?. No delay greater than 
3 10-1! sec at this point was introduced by magnetization. 

The long period in distrene reported by Bell and Graham (1952) was 
investigated further by the comparison method using aluminium as 
standard. 

482 
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At this stage, a new ?2Na source was prepared for these experiments, 
by boiling up turnings from a deuteron bombarded magnesium target, 
kindly supplied by Professor Burcham of Birmingham University. As 
the magnesium was oxidized the sodium went into solution as NaOH and 


Fig. 2 


No. of Coincidences. 


Delay x 10-*s. 
Coincidence resolution curve of coincidences between the 1-28 Mev y-ray of 
72Na in C.1 and 511 kev annihilation y-rays in C.2. 


© positrons annihilating in lead. 
x positrons annihilating in iron. 


was readily separated from the MgOH by virtue of its vastly greater 
solubility. The solution was deposited on 0-0002 inch aluminium foil, and 
covered with a foil of similar thickness. The total weight including 
backing was 4 mg/cm?. 

The discriminator D.1 (see fig. 1) was set so as to exclude 511 kev 
y-radiations and the discriminator D.2 was set so as to exclude the 
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23H g 280 kev y-radiations. This helps to reduce spurious coincidences 
due to scattering. Care was taken with the geometry for this reason. 
A resolution curve of positrons in distrene and a comparison enrve using 
aluminium next the source were taken at room temperature (15°C). The 
curves obtained are shown in fig. 3. The detailed shape of the resolution 
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Delay x 10-°s. 
Coincidence resolution curves of coincidences between the 1-28 Mev y-ray of 
22Na and 511 kev annihilation y-rays in C.2. 


positrons annihilating in distrene. 
© positrons annihilating in aluminium. 


curve for distrene, the large difference between the lifetime measured from 
the tail and from the centroid shift, showed that the long lifetime was 
associated with about one-third of the positrons. Assuming that the 
shorter period was comparable with that of aluminium the analysis of 
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Bay (1951) using second moments was applied. This gave the second 
mean life—2-5-++0-3 x 10-® see for a fraction 0-27--0-03 of the positrons. 
This is in good agreement with the period as found from the tail. These 
results agree with the recent work published by Bell and Graham (1953). 
The general features of the curves also correspond to those of the above 
authors. 

To enable variation of lifetime with temperature to be studied without 
affecting the apparatus, the scintillators were surrounded by a jacket 
through which cold air was blown. The specimens were electrically 
heated. 

The variation of lifetime with temperature was examined at a number 
of temperatures between room temperature (15°C) and 130°c when distrene 
begins to soften. An increase of about 10% was found for the 115°c 
increase in temperature. 

The internal structural parameters of plastics are thought to vary with 
method of preparation and treatment. Samples machined from rod, cut 
from 0-1 inch sheet, made up from layers of 0-002 foil were examined. No 
difference was found between lifetimes in any of these. 

No residual effects of heating, even to temperatures at which the distrene 
became soft, were found. 

Plastic sulphur was examined in the same way as distrene. The 
samples, prepared by distillation on to the cold aluminium plates used in 
the experiment, were 0-1 inch thick and of a transparent straw colour. 
To avoid physical change the experiments were carried out immediately 
after preparation. With the resolution used no evidence indicating a 
second period exceeding 8 x 10-19 sec was found. Bell and Graham have 
reported some evidence for a partial period of about 4x 10~?° see. 


§ 5. SEARCH FOR THREE-QUANTUM RADIATIONS FROM DISTRENE 


On the hypothesis that the long period in certain materials is associated 
with three-quantum annihilation from the triplet state of positronium a 
search for this was made by two methods. In these, sodium iodide 
counters were used, with coincidence units of 0-5 see resolving time. 

The first method used three counters at 120°, arranged to minimize 
scattering, and a ®Cu source 0-0001 inch thick surrounded by distrene. 
The discriminators let pass only radiations above 280 kev. After 
subtracting randoms the number of coincidences was 


Counters Coplanar 25-+-5, 
NON 13+3:5. 
A check experiment was made using aluminium in place of distrene. It 
is seen that little three-quantum annihilation occurs. 
The second method measured the y-ray coincidences at 180° using two 
counters 20 cm apart and the thin Na source at the mid point of their 
axis. This geometry favours detection of two-quantum as against 
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three-quantum annihilations by 100:1. The number of coincidences 
was measured with distrene and aluminium as absorbers and results 
indicated that less than 2°, of the annihilations are three-quantum. 


§ 6. LireTmmMs or Postrrons In ALUMINIUM 


To measure the lifetime of positrons in metals it is necessary to produce a 
resolution curve for prompt events of similar energy. From the spins and 
parities of the two y-emitting states in °°Co (cf. for instance Metzger and 
Deutsch 1952) the delay in emission for the second y-ray can be expected 
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to be about 10-!2 sec (cf. Moon 1949). In the case of *?Na the 1:28 Mev 
y-ray is thought to be prompt. An electric quadrupole transition would 
again give a lifetime of about 10-”sec. A transition of higher order 
would have a lifetime readily detectable here. The experiments should 
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thus give an absolute measurement of the positron lifetime within the 
errors quoted below. 

Those parts of the ®°Co y-spectrum and the 1-28 Mey *Na y-spectrum 
which lay above the bias of discriminator D.1 were examined with a 
kicksorter and found to lie close to one another. 

Instead of counting all pulses from C.2 exceeding D.2 bias, only a band 
of pulses was accepted. The band was positioned at a region where the 
annihilation radiation spectrum and the ®°Co y-ray spectrum were similar 
in shape. Thus in each counter the accepted spectrum of ??Na and ®°Co 
was the same and hence their resolution curves could be compared. 

In the first experiments the band was chosen by discriminator D.2 and 
another in parallel with it connected by an anticoincidence unit. Later 
the full apparatus shown in fig. 1 was used. The pulses fed to the kick- 
sorter were about 10-7 sec wide. Delay of 1-25 usec was inserted as 
shown using Z.14M delay cable so that these might coincide with the 
corresponding gateing pulses from 3-C. 

Using the thin ??Na source described above, surrounded by aluminium, 
a resolution curve was taken. The experiment was repeated using a ®°Co 
source of similar strength, giving a prompt resolution curve for comparison 
(see fig. 4). Several runs, sometimes in succession, were taken in this 
manner using both arrangements described above. Good agreement was 
obtained between their results. Each pair of curves was obtained in about 
two hours, with counting rates of a few hundred a minute at peak. Errors 
were mainly due to possible energy drifts of side channels. The result 
obtained was (1-6-L0-6 x 10~1° sec). This is in good agreement with the 
value obtained by Bell and Graham (1953) using a different method. 

A possible development of this technique may be mentioned, particularly 
useful for investigating delayed y-cascades. By selecting only a suitable 
band of pulses from counter 1 also, prompt resolution curves from say the 
y-y coincidences of annihilation radiation could be obtained. These 
could be compared with a wide range of suspected delayed resolution 
curves. 


§ 7. Discussion 


The measured lifetime in aluminium, and by the comparison experi- 
ments, in all metals so far examined, is much shorter than that expected 
from simple annihilation of positrons. It would appear that the cross 
section for positron capture may be considerably greater than that for 
annihilation and that most annihilation comes from bound states. These 
considerations are strengthened by the results of De Benedetti et al. (1950) 
and Du Mond, Lind and Watson, which suggest that only the valence 
electrons take part in the annihilation process. Again, alignment of the 
d-shell electrons of iron does not seem important. 

The experiments on distrene show that if the second period is due to 
triplet positronium, much of it is rapidly converted to the singlet state. 
If the long lived fraction were to find its way into sheltered states such as 
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“vacancies in the matrix, the interaction might be related to the spacing 
and so decrease with temperature. 
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SUMMARY 

Experiments are described in which measurements of pressure gradient 
were taken during the flow of mercury in a rectangular channel of 15: 1 
aspect ratio, in the presence of a transverse magnetic field. It was 
possible to suppress turbulence even at the highest obtainable Reynolds 
Number (R=10°). 

For high values of R, the friction coefficient C; appears to be a function 
of M/Rand fp. It is suggested that when R is large enough C; might be a 
function of the single group M/RB'2. This is comparable with the simple 
relationship proposed by Hartmann and Lazarus (1937). 

Lundquist (1952) has suggested that the transition from laminar to 
turbulent flow occurs at a fixed value of R/M, and that 6 does not enter 
into this criterion. The experiments suggest that the critical condition 
for the turbulent to laminar transition is R/M ~ 900, although since only 
one liquid was used, it was not possible to determine the influence of f. 


§1. THe Work oF HARTMANN AND LAZARUS 


THE equations of motion for fully developed laminar flow between parallel 
plates, of a viscous, electrically conducting fluid, in the presence of a 
transverse magnetic field have been solved by Hartmann (1937). Accord- 
ing to this theory the magnetic field flattens the velocity profile and 
increases the rate of shear at the walls, thus increasing the pressure 
gradient necessary to obtain the same rate of flow. 

Hartmann and Lazarus (1937) carried out a series of experiments on the 
flow of mercury in rectangular and circular tubes with an imposed trans- 
verse magnetic field. They produced a series of curves showing pressure 
gradient as a function of field strength, each curve representing a constant 
volumetric flow rate. With flow which was laminar in the absence of a 
field, the effect of the field was to increase the pressure gradient in the 
manner predicted by Hartmann. When the flow was turbulent in the 
absence of a field, increase of field strength produced a decrease in pressure 
gradient up to a certain point, after which it increased again, finally 
obeying Hartmann’s theory. 


* Communicated by the Author. 
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The experimenters assumed that the phenomena observed in the 
turbulent regime were due to this ‘ Hartmann Effect ’ combined with the 
progressive damping of turbulence by the field. In order to separate these 
effects, they assumed that turbulent damping was proportional to H?, 
where H is the field strength. On this basis, and after making other 
assumptions, they deduced that the decrease in pressure gradient due to 
damping of turbulence, 4),, was proportional to H2/d, where d is the tube 
diameter, and was independent of the mean velocity u. They sought to 
confirm this by a dimensionless analysis assuming that the decrease in 
pressure gradient was a function of H, d and w. 

When conditions do not change with time, and when there are no free 
surfaces, magneto-hydrodynamic channel flow of a liquid metal depends 
upon the seven quantities H,d,u,p,y, 0,4, where d is a characteristic 
length (taken as the semi-width in the case of a high aspect ratio channel), 
p.7, o, » are respectively the density, viscosity, electrical conductivity and 
permeability of the fluid. The following dimensionless groups are parti- 
cularly useful in channel flow experiments : 


4oud 47 
jpeg m=und|(2), (jae 
7 7 p 


f& is the usual Reynolds Number, and f is a fluid property whose signifi- 
cance for the amplification of a small magnetic disturbance has been 
pointed out by Batchelor (1950). £ is the ratio of two diffusivities, and 
in this sense is analogous to a Prandtl Number. 

The dimensional analysis carried out by Hartmann and Lazarus tacitly 
assumes that 4, is independent of 7 and o. When the flow is highly 
turbulent (their flow was not highly turbulent) it is plausible by analogy 
with ordinary hydrodynamic flow to assume independence of 7, but they 
had no evidence to assume independence of o. 


§2. EXPERIMENTAL METHODS AND RESULTS 


2.1. Hxperimental Equipment 


The experiment to be described was designed to study the relationship 
between the pressure gradient, mean flow velocity and field strength, in 
the turbulent flow of mercury in a channel in the presence of a transverse 
magnetic field. The previous work of Hartmann and Lazarus was carried 
out with Reynolds Numbers below 5500, and so is hardly out of the 
transition range. The present apparatus was designed to give Reynolds 
Numbers up to 10°, and at the same time a much higher value of M than 
was possible in the apparatus of Hartmann and Lazarus. 

The channel was constructed from non-magnetic stainless steel plates, 
insulated internally with a coating of ‘Araldite’. Its length was 
36in. and the internal cross section measured 4in.x0-266in. The 
width (0-266 in.) was checked at many points with an internal micrometer, 
the departure from the mean being less than 0-003 in, Over the test 
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section the width was constant to within 0-:001in. Twenty-one static 
pressure holes were carefully drilled and broached in the top of the channel. 

Pressures were measured on a multi-tube manometer. Since a gravity 
feed was used, the mean static pressure in the circuit varied during each test. 
For this reason the upper ends of the manometers were connected via a 
manifold to a small closed tank, the entire volume above the mercury 
being filled with air-free distilled water. 

The volumetric flow rate was measured with a magnetic flowmeter, 
consisting of a section of insulated tube placed in a transverse magnetic 
field. Two electrodes were placed in the tube at opposite ends of a diameter, 
this diameter being perpendicular to the directions of flow and field. The 
flowmeter output was measured on a potentiometer, and the device was 
calibrated by measuring the throughput in a given time. The accuracy of 
calibration is believed to be +14%%. 


Fig. 1 
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The fluid circuit is shown schematically in fig. 1. Freshly distilled 
mercury was contained in four oxygen cylinders placed horizontally about 
20 ft. above ground level. On opening valve (a) flow was established 
through the flowmeter and channel into the sump tanks. At the end of 
each test, the sump was pressurized and mercury returned to the header 
tanks via the valve (b). The 15° inclination of the channel, in conjunc- 
tion with the stand pipe, ensured that at the end of each test the channel 
and manometers remained full of mercury. The whole circuit was 
designed so as not to trap any air, since this would eventually displace 
mercury from one of the manometers. Great difficulty was experienced 
in removing all traces of air. 
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The transverse field was provided by a calibrated electromagnet having 
pole faces 12 in. x5 in. and a gap of lin. The flux distribution in the gap 
was checked and found to be uniform within gale ay 

The last 12 in. of the channel were located in the magnetic gap, leaving 
a 24 in. entry length for the full development of turbulent flow. 


2.2. Hxperimental Results 


During the tests it was observed that the pressure gradient remained 
constant over the central six inches of the field, and increased rapidly at 
the edges. This was presumably due to the edge effect currents mentioned 


Fig. 2 
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by Hartmann and Lazarus. Since only a limited time was available for 
each test, it was considered more profitable to concentrate on the central 
region of constant pressure gradient. Manometers 5 to 11 were read in 
semi-random order, care being taken to obtain repeated readings of at least 
numbers 5 and 11. Typical pressure gradients are plotted in fig. 2. Where 
only one point is given opposite manometers 5 and 11, the repeated readings 


coincided. 
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The table shows the complete results in dimensionless form. Cy is the 
usual friction coefficient, defined by 7y=C,pu?/2, where 7) is the wall 
shear stress. 


LO-h M 10°C; 10-5F M 108Cy 
1-21 0 4-27 0-88 120-5 11-2 
1-22 22-0 4-32 0-87 128-8 12-0 
1-21 31-0 4-48 0-80 0 4-81 
1-21 60:3 531 0-80 11-5 4-65 
1-21 66:8 5:58 0-81 24-2 4-87 
1-22 86:5 6-39 0-80 35°9 5-31 
1-21 99-8 7-02 0-80 46-2 5-74 
1-2] 112-5 7:67 0-80 51:9 6-03 
1-21 125-5 8:33 0-80 80:3 8-11 
1-21 129-0 8-80 0-81 115-0 11-73 
1-01 0 4-58 0-70 0 4-99 
1-02 14:3 4-41 0-70 20°3 4-90 
1-02 32:6 4-72 0-70 29-0 5:00 
1-02 40-7 4-98 0-705 42-4 5°82 
1-01 49-2 517 0-705 48-4 6-29 
1-01 72-2 6-41 0-70 60-0 7:24 
1-00 105°8 8-51 0-70 110-5 13-81 
1-00 118-5 9:36 0-70 124-7 14-32 
PAVL 130°5 10-32 0-30 0 5-95 
0:88 0 4-67 0-30 5-2 5-81 
0-88 0 4-63 0-30 10-4 5°86 
0-88 23:0 4-72 0-30 18-0 6-01 
0-88 30-0 4:87 0°30 26°0 7°05 
0-88 42-1 5°39 0-30 63°1 17-1 
0-87 75:7 7:38 0-10 0 7:73 
0-87 99-3 7-49 0-10 5-4 5-67 
0-87 105-8 9-85 0-10 16°4 9-27 

0-10 46-0 37-0 


$3. Discussion 


Figure 3 shows a graph of CO; against M/R and shows several interesting 
features. Apart from the two low Reynolds Number curves (which are of 
lower accuracy owing to the small pressure gradient) the experimental 
points lie quite well on a single curve, except at low values of M/R. This 
curve approaches the laminar ‘ Hartmann’ line slowly, which seems to 
indicate that turbulence is not suddenly suppressed, or that if it is, it 
persists so long that its effect on C; is swamped by the Hartmann effect. 
Within the experimental scatter, the flow appears to be laminar whenever 
R<900M. It is interesting that laminar flow can be obtained in this 
manner even when R= 105. 

Figure 3 shows that at the higher Reynolds Numbers, the trend is for 
C, to be a function of M/R, and that the range of validity of this relation- 
ship widens with increasing R. One may speculate that at high enough 
values of R, C, is a function of M/R and B only. Now if at high values of 
hk, the effects of viscosity can be neglected by comparison with the other 
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forces present, then dimensional reasoning shows that C, must be a 
function of the two groups (M/RB'!2) and (RB). If this is to be a function 
of (M, /R) and £ then this function must be of the single group (M/Rp'!). 
It is interesting to note that this quantity does not contain dissipative 
terms and is, in fact, equal to the square-root of the ratio of magnetic to 
kinetic energy per unit volume. This speculation is a generalization of 
the relationship intuitively assumed by Hartmann and Lazarus. They 
assumed a simple form for the function, and moreover did not restrict it 


Fig. 3 
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to high values of #, although their attempts to separate the Hartmann 
effect were equivalent to removing direct viscous action. They had no 
justification for omitting o from their dimensional analysis. 

At first sight there appears to be a paradox since (M/Rf'!) does not 
contain c, and so the speculations might equally apply to a non-conducting 
fluid. It must be remembered, however, that JM and F# can separately 
enter into the conditions under which the speculations apply. These 
conditions remain as yet undetermined. 

In a recent paper (Lundquist 1952) it is suggested that transition from 
laminar flow occurs at a certain critical value of the boundary layer 
Reynolds Number R’ ,,i,(R’=4puL'/n, where L’ is a dimension character- 
izing the boundary layer). In the present notation L’=d/M and thus 
B’ i= (R/M) exit, and is independent of f. It appears from fig. 3 that 
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transition the other way (i.e. from turbulent to laminar flow) takes place 
at approximately R/M=900, although since only one fluid was used, it is 
not possible to assert or deny independence of f. 

The essence of Lundquist’s argument is that magneto-hydrodynamic 
forces are important only insofar as they determine the thickness and 
velocity distribution in the boundary layer, and that apart from this, they 
do not enter into the stability criterion. The parameter f, therefore, does 
not affect stability. It is very desirable that similar experiments be 
carried out using other fluids, and so varying f. 

The present writer has attempted to separate the Hartmann effect from 
the turbulent damping effect in a more refined way than that of Hartmann 
and Lazarus. The essence of the method is to apply Hartmann’s equations 
to a ‘magneto-hydrodynamic boundary layer’, i.e. a layer near the 
wall in which only electro-magnetic and viscous forces are important. 
A relationship follows between the shear at the two boundaries of this layer 
in terms of M, R and the thickness of the layer. Thus knowing the 
velocity gradient at the wall, from the pressure gradient, the velocity 
gradient at the inner boundary of the layer can be estimated, and this can 
be used to measure the dissipation in the turbulent core. The results 
obtained, however, are very sensitive to the choice of boundary layer 
thickness, which at present is arbitrary. Until more experimental 
evidence is available it is not desirable to publish these results. 
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SUMMARY 


The idea of conduction by electrons in the energy bands of an impurity 
system isexamined. For a lattice of hydrogen-like impurities the edges of 
the Is, 2s, and 2p bands are calculated. Effective masses for electrons 
in these bands are introduced. With their help measurements of the 
conductivity and the Hall constant are interpreted. Depending on the 
density of the impurities and on the temperature, conduction takes place 
predominantly either in the band of the medium, or in excited impurity 
states, or in the 1s impurity band. In this last case the activation energy 
vanishes. 

This interpretation implies a description of the electronic state of the 
impurities with band wave functions and indicates the validity of such a 
description up to rather large values of the lattice constant. This can, 
however, be reconciled with the insulating properties of the oxides of 
transition metals. 


§ 1. INTRODUCTION 


In recent years various authors have indicated the importance of 
conduction by charges in impurity states. Busch and Labhart (1946) 
explained their measurements of the conductivity and the Hall constant of 
SiC with a model which attributes a mobility to the electrons in impurity 
levels. Similarly Hung (1950) postulated conduction in impurity bands in 
view of Hall effect measurements on Ge at low temperatures. Erginsoy 
(1950, 1952) suggested an investigation of impurity bands and calculated 
the lower edges of the s bands. Earlier the work of Pearson and Bardeen 
on Si (1949) had shown that the activation energy of the charges depends 
on the impurity content, and these authors had assumed that this results 
from an interaction between the free charges and the ionized impurities. 
However, this argument was criticized by Castellan and Seitz (1951), who 
investigated various alternative theories without finding a satisfactory 
explanation. In particular, for a density of impurities such that the 
samples show a metallic characteristic the interaction energy between the 
impurity centres in their ground state turned out to be small. Hence the 
Heitler-London scheme seemed to be the appropriate description of the 
impurity system, whereas the Bloch scheme which could describe a metallic 
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behaviour was considered as representing an excited state. The possibility 
had to be faced that fundamental changes in the theory would be necessary. 

We take up the idea of conduction in impurity bands without justifying 
the application of the band scheme. First we investigate the upper and 
lower edges of the 1s, 2s, and 2p bands, and then define an effective mass 
for an electron in each of these bands. A comparison with the effective 
mass for the conduction band of the medium suggests an interpretation of 
conductivity and Hall effect measurements, in particular the observed 
decrease of the activation energy. The agreement with experfment is 
significant for solid state theory, since it indicates the applicability of the 
band scheme up to rather large values of the lattice constant. 


§2. THe ENERGY Banps 


As a model for a semiconductor we take a medium, characterized by the 
effective mass m* of a conduction electron, and by the dielectric constant x, 
sufficiently large, in order that the hydrogen-like picture of a single 
impurity is justified (Baltensperger 1951). We assume that the impurities 
are arranged in a regular close-packed lattice. Applying the Wigner— 
Seitz method atomic polyhedra are set up, which we approximate by 
spheres of equal volume. Their radius 7, is then related to the density 
of the impurities V by 

(47/3) r 8=1/N. . oo OS 


On the average an electron in a sphere will be subject only to the ionic 
potential V=—e/xr. The Schrédinger equation 


h? e 
Ime 4b (= +«) Y=0 2 6 <2 ets, ee 
has the general regular solution 

p=, (7r)X Yi, »(% ?) Me. Oe ee 
FR, (r)=exp (—1/nxa) x (2r/nxa)'x F(l+-1—n, 214-2, 2r/nea),. (4) 
where Fey alate + OS + me | (5) 

is the confluent hypergeometric function, and 
a=h?lm*e? eee SoA ry a 

The energy is given by e=— UA & = 
Wee * pe? + 0 ie ee 


where n is to be determined by an appropriate boundary condition. 
This boundary condition is obtained from the requirement that the 
wave functions have Bloch’s form (Slater 1934) 


b=exp (27ikx) x u(x) 
u(x-+-t)=u(x) oleae ee he, feb, ae ec 
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where t is a primitive translation of the lattice. We only consider special 
points of a reduced Brillouin zone, namely the centre, k=0, to which 
corresponds a periodic wave function, and the corners, given by 2(kt)=1, 
the wave function thus being antiperiodic. With the spherical approxi- 
mation for an atomic cell, J is a good quantum number. If a wave 
function is even about the nucleus (even J), the corresponding boundary 
condition is Or) 

AUS) eae led eee ae em Pn 6) 


or T=? 

ame | 

for k=0 and Win) = Ua eS ee hye eee (1) 
for 2(kt)=1. For odd wave functions k=0 belongs to (10) and 2(kt)=1 


to (9). We shall consider the eigenvalues belonging to (9) and (10) to be 
respectively the bottom and the top edges of the bands. 


Fig. 1 
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Using (3) and (4) the conditions read 
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Equation (10a) has been used by Michels, De Boer and Bijl (1937), 
Sommerfeld and Welker (1938), and by De Groot and Ten Seidam (1946) 
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to determine n, and therefore «, as a function of r./Ka and J. Further, 
eqn. (9 a) has been discussed by Erginsoy (1952) with =0, and n approxi- 
mately equal to an integer. Similar methods can be used for J=1 and 
some further values can be obtained from existing tables of F(a, Y x) 
(Webb and Airey 1918, Jahnke and Emde 1945). Forno the conditions 
(214.2) x (89/1) 2? X For4s([8r [xa] =Fo([8r/ea}*?) . . . (98) 

J oa ((8rq{ea)/2)=0 . 2 1 


give the values of r,/«a for e=0. Theedges thus obtained are shown in fig. 1. 


§ 3. DEFINITION OF EFFECTIVE MASSES 


Instead of discussing the band widths de, we introduce quantities m* 
defined by 
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where & is the largest wave-number for a spherical Brillouin zone and 
satisfies 

(4n/3)S=N. 2. 2 
From (11), (12), (1) and (6), we get 


m* ( 4 \18 aes r\2 de 
aor rs 32 ~ pen ve Mde=0-1710 (=) < eka’ ; (13) 
For r>0, m*/m' was found by solving (2) with (9) and (10), neglecting the 
potential energy. Introducing this as a perturbation, we obtain the slope 
of m*/m' forr,>0. In fig. 2, m*/m' is plotted against r,/«a for the 1s, 2s, 
and 2p bands. 
While m* characterizes the dynamical property of an electron moving 
through the medium, either free or bound to an impurity, m" will be 
considered as the effective mass of an electron whose energy lies in an 
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impurity band, for its movement through the impurity lattice. The 
definition (11), (12) of a constant m* for a whole band is of course rough, 
since it implies a parabolic dependence of the energy on the wave-number 
for an entire Brillouin zone. 

The mobility of an electron is finite because of deviations from a perfectly 
regular lattice. For an electron in the conduction band scattered by 
thermal vibrations, the mobility is proportional to m*-*/2 (Seitz 1948), 
while scattering due to impurity ions according to the formula of Conwell 
and Weisskopf (1946, 1950) gives a dependence like m*-1/2. The conduction 
mechanism of electrons in impurity states is obscured by the fact that the 
impurities in reality do not form a lattice. On the one hand irregularities 
in the arrangement of the impurities will lead to scattering. On the other 
hand the disorder will increase the average interaction between the 
impurities, which is a minimum for a regular arrangement, so that the 
present model overestimates the effective mass. This is confirmed by 
a calculation of Aigrin and Jancovici (private communication) who have 
calculated the energy levels for a disordered arrangement of the impurities. 
The exact relation between the mobility of an electron in an impurity 
state and the effective mass m' is not known. For a qualitative inter- 
pretation of the measurements it is sufficient to assume that m*/m' 
(fig. 2) is of the order of magnitude of the ratio of the mobility belonging 
to the impurity band to the mobility of the band of the medium. 


§ 4. INTERPRETATION OF MEASUREMENTS 


The temperature determines the distribution of the charges amongst 
the bands. If there are NV, charges in the ith band having a mobility 
4;, the conductivity is 


Oto Ne feps we ee es (14) 

and the Hall coefficient (in prachical units) 
R=F (Xe; Nu ?)x(&e; NOM ee toes ee (15) 
where ¢,—-e if the cerees are holes, ‘sii e,=—e if they are electrons, and 


where F is a constant of order 1. If only one band, the ith say, contributes 
to the conduction, the Hall coefficient determines the number of charge 


carriers : R=F(e,N;)7 if Nj pj=9 for (Ean So ep MP ao (16) 


We consider the measurements of Pearson and Bardeen (1949) (see table 
1). For sample 1 the mobilities of the impurity bands are negligible, hence 
the measured activation energy is the energy required to raise a charge into 
the band of the medium. In samples 2, 3, 4, A, B, activation into an 
excited impurity band is sufficient to provide conduction. Sample 5 
with r,/ka=5-2 is intermediate. The mobility in the 1s band is small but 
not negligible. Thus at high temperatures the resistivity of sample 
5 behaves similarly to that of sample 4, while at low teniperatures the 
conduction in the Is band becomes predominant, and then the measured 
activation energy vanishes. Samples 6, 7, 8, C, D, have the metallic 
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behaviour corresponding to conduction in the Is band. The Hall 
coefficient shows that in these cases there is one charge carrier per impurity 
atom. The mobilities are largest for the samples 2, 3, 4, A, B, for which 
the excited impurity bands are important. However, it must be 
remembered that when different bands contribute to the conduction the 
determination of the numbers of charge carriers, the mobilities, or the 
mean free paths becomes ambiguous. ; 
Bands belonging to large quantum numbers n have not been considered. 
Since the extension of hydrogenic wavefunctions is essentially proportional 
to n, we expect by extrapolation that for intermediate temperatures these 
bands are important in a range, 6n>r,/ka>2n. In fact, for low densities, 
r|ka>6n, the mobilities will be negligible, while for high densities, 
7,|Ka<2n, the bands are shifted to such high energies that only few 
electrons are activated into them. , 


Table 1. Interpretation of the Measurements of Pearson and Bardeen 
(1949) on Silicon Alloys 


eT ee ee 


Carriers/cc. Conductivity 
Sample Solute in saturation 1's/KQ mainly due to 
range 

1 Boron 8-0'< 10%* 45 Silicon band 

2 Boron 6-0 x 101" 10-7 Excited impurity 

3 Boron F3x 10% 8:3 bands 

4 Boron 2°5 < 1048 6-6 (28, 2p, 3s, 3p, 3d, . .) 

5 Boron 5-3 x 1018 5:2 1s band for low tem- 

perature 

6 Boron 1-4 x 1019 3-7 

7 Boron Tax Lote 1-8 ls band 

8 Boron 4-8 x 10° 1-15 

A Phosphorus 1-05 x 10?” 13 Excited impurity 

B Phosphorus 120261048 5-6 bands 

C Phosphorus Lie OFe 2-4 ab 

D | Phosphorus | 2-6 1020 0-05. || faBsPane 


rs/ka according to (1) and (6) with «=13, and 
m*—m for the silicon—boron alloys, 
m*=0-67 m for the silicon-phosphorus alloys. 


§ 5. SIGNIFICANCE FOR Souip Starr THEORY 


Although the disordered arrangement of the impurities makes quantita- 
tive predictions difficult, experiments with semiconductors may indicate 
the qualitative behaviour of a system of atoms for a large range of 
densities. In the usual theoretical method the electron state is approxi- 
mated by a determinant of one-electron functions. Tn the Heitler—London 
scheme these are localized, and, for a description .of the ground state of a 
system with NV ions and N electrons, N different spatial functions will be 
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occupied. In the Bloch scheme the functions extend over the whole 
lattice, and in the ground state N/2 functions belonging to the lowest 
energies will be occupied twice. 

In fig. 3, as an illustration, the total energy of the system according tothe 
two schemes is plotted against r,/«a. For the Heitler-London scheme the 
curve corresponds to a formula which was derived by Slater (1930) 
disregarding the normalization difficulties that arise from non-orthogonal 
wavefunctions. However, for a system containing a large number of 
electrons, this invalidates the quantitative results of such a calculation 
(Slater 1930, 1953, Loéwdin 1950). As to the Bloch scheme, for high 
densities the curve represents the results of the calculations on metallic 
hydrogen by Kronig, de Boer, and Korringa (1946), while for large inter- 


Fig. 3 
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total energy for the Bloch scheme. 
total energy for the Heitler-London scheme. 


atomic distances it is based on the considerations by Seitz (1940). The 
occupation of the N/2 functions with lowest energies is a good approxi- 
mation for wide bands only, as was pointed out by Slater (1953), whereas 
for narrow bands the wavefunctions should be chosen to be a sum of 
determinants corresponding to all different configurations. A calculation 
of this type may yield a much lower energy. At present no reliable 
quantitative treatment is available for a lattice of atoms with weak 
interaction. 

The two types of description are qualitatively different (Schubin and 
Wonsowski 1934). At low temperatures the Heitler—London scheme does 
not allow any conduction. This is true provided the numbers of electrons 
and atoms are exactly equal; otherwise the difference between these 
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numbers gives the number of free charges. The Bloch scheme with 
incompletely filled bands describes a metallic conductor. ee. 

Mott (1949) has pointed out that the state of a substance will either 
correspond to one or the other scheme, and not to some intermediate 
description. Oxides of transition metals, such as NiO, are insulators in 
the pure state, although they have partially filled energy bands; for 
these the Heitler-London scheme must apply. If the lattice constant of a 
metallic substance is increased, at a critical value a phase transition will 
occur, corresponding to a change from the Bloch to the Heitler—London 
description. For impurities in a disordered arrangement this transition 
may not be sharp (Mott 1951). 

We have seen that the idea of conduction in impurity bands explains 
the dependence of the conductivity and the Hall coefficient on impurity 
content and temperature. It implies that the impurity states are des- 
cribed by band wave functions, whenever they are found to contribute to 
the conduction. It should be mentioned that a finite conductivity at 
low temperatures alone could arise from some secondary effect ; since 
a sample will contain both donors and acceptors, and since one impurity 
type will trap charges of the other, even the Heitler—London description 
permits conduction, though for this the number of free charges is given by 
the density of the secondary impurities. However, the measurements of 
Pearson and Bardeen (1949), and, more especially, those of Hung and 
Gliessmann (1950) show that at low temperatures, at which conduction in 
excited states becomes negligible, the Hall coefficient is constant and 
approximately equal to the Hall coefficient at saturation temperatures, at 
which all electrons are in the conduction band of the medium. 

This observation demands a description of the ground state of the 
impurity system in which all electrons are free to move (Hung 1950). It 
excludes an explanation of the low temperature conductivity by a 
secondary effect. It also indicates, that if the phase transition to the 
Heitler-London scheme is not sharp, it is not yet noticeable in the 
samples, for which the constant Hall coefficient at low temperature 
was measured. Sample 5, table 1 then gives a lower limit for the critical 
radius of the atomic sphere for which the phase transition occurs 

r[Ka>5. 2 «+. se) ee 
The value r,/ka=2-5 of fig. 3 shows that the underlying theory is 
inadequate, even considering that the disordered arrangement of the 
Impurities is a source of error. 

The question arises, whether the condition (17) can be reconciled with 
the insulating properties of the oxides of transition metals, which demand a 
Heitler-London description for the 3d states of the metal ions (Mott 1949). 
The above discussion was based on the parameter r ,/Ka, ie. the radius of 
an atomic sphere measured in units which are atomic units divided by the 
effective number of charges, I/k. For a rough argument we define a 
similar parameter for metal oxides, such that for equal values of the para- 
meters we may expect an interaction of comparable relative importance. 
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NiO, for example, has a simple cubic structure, in which the Ni** ions 
alone form a cubic face-centred lattice. The atomic volume is 18-2 AS 
so that r,=1-63 4. We express r, in units NOo|Zegr, Where a)=0-529 A 
is the Bohr radius, and n=3 accounts for the increase in size of the 
wavefunctions with higher quantum numbers. The effective number of 
charges Z. for one of the 8 electrons in the 3d shell of a Nit+ ion will be 
larger than 3, since each of the remaining 7 electrons in the shell screens 
the nucleus only for a fraction of the time. If this fraction is 4, then 
Zog=3+4X7=6-5 and | ; 
1 Loe  163X 65 

Ny 3X0-B29 
The value of the parameter is larger than the minimum required by 
(17). It must be remembered, however, that this argument ignores the 
presence of the O-~ ions, and also that for very large values of the para- 
meter even the Bloch scheme leads to an unobservably small conductivity. 

Thus the evidence from semiconductors and from the oxides of transition 
metals are consistent with the view that a description of a solid must be 
either of the Bloch or the Heitler—-London type but the phase transition 
between the two seems to take place for larger values of the lattice constant 
than expected according to the present simple theory. 


(18) 


I am grateful to Professor L. Rosenfeld for his kind support and 
stimulating discussions. 
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ABSTRACT 


The band theory of metals is applied to analyse the thermoelectric 
power of ‘simple’ metals at high temperatures. The conventional 
analysis is first carried through for entirely free electrons and then 
extended to include variation in number of free electrons per atom, 
departure from a quadratic dependence of energy on wave number, and 
non-spherical Fermi surface. Consideration of the data on the alkali and 
noble metals shows potassium to present a special problem in interpre- 
tation. 


§ 1. INTRODUCTION 


On the basis of the electron theory of metals, the variation with tempera- 
ture of the electrical and thermal conductivity can be discussed, and the 
theory has been used to study particularly the alkali metals, which are 
often regarded as more or less ideal. In the case of electrical resistance, 
the features of the dependence on temperature are quite well predicted, 
but the agreement between theory and experiment is less satisfactory 
for thermal conductivity. The actual magnitude of the conductivities 
is more difficult to analyse without an explicit knowledge of the potential 
energy of an electron in the metallic lattice. 

The phenomenon of thermoelectricity presents fresh problems. At 
low temperatures, such that impurity scattering becomes significant, a 
rather detailed knowledge of the mechanism of scattering of electrons 
by impurity atoms appears to be called for to explain the varied anomalies 
in thermoelectric behaviour in this region (e.g., Borelius 1935, MacDonald 
and Pearson 1953). Such theories of impurity scattering as are available 
at present are inadequate to account for the experimental results. 

In this paper we shall confine our attention to the region of ‘ high’ 
temperatures, that is, when 7’, the temperature of the metal, is at least 
of the order of the Debye characteristic temperature 9, but is below the 
melting point. 

This region of temperature is then classical in nature, a relaxation 
time always exists, and impurity scattering is in general negligible. 
Furthermore, the Einstein model to describe the heat motion of the 
lattice should be adequate. Rather definite theoretical conclusions may 
then be derived on the basis of this model, and compared with experi- 
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mental results to throw further light on the basic assumptions of the 
theory. 


§ 2. EXPRESSION FOR THERMOELECTRIC POWER 


The absolute thermoelectric power S of a metal is given by (for a fuller 
account see Wilson 1936, Mott and Jones 1936) 


s= (3-2) Shee ie Pe Aer 


where ¢ is the Fermi energy at the temperature 7’, ¢ is the charge on the 
electron, and K,,, is given by 


z Te 0H \ 207, 
K.=— go | | 2" (ae) ap r(kidk,dkydk,. . . . (2) 


The integration in (2) extends effectively throughout the volume enclosed 
by the Fermi surface, and /, is the equilibrium distribution function. ZF is 
the energy of the electron in the state defined by the wave number k of the 
electron, k,, is the component of k in the x-direction in which the external 
field and the temperature gradient exist, and 7(k) is the relaxation time. 
It is not assumed here that 7(k) is a function of | k| only as is done in 
the case of electrical and thermal conductivity. Equation (1) should, 
therefore, be valid for any relation between H and k (cf. Mott and Jones, 
esl): 
Substitution of (2) in (1) leads to 


w2k?T ( o 7(k) OE\2 


which is the expression given by Wilson (p. 176). The integration is 
over the surface of constant energy #; E, is the Fermi energy corres- 
ponding to 7’=0. 


Writing now ev | OEN? __7z(k) 
CL y= i iiaradiely (eb ee ewe ent 4) 
in eqn. (3), we get finally 
__ kT [0 log o(£) (5) 
ely \Olog Ff |poz. 


which is the expression given by Mott and Jones (p. 310). 
For perfectly free electrons, the surfaces of constant energy are spheres. 
If we then take 7(k) as a function of k only, eqn. (4) can be written in 


the form o(£)=const. | k |°7, aoe Pans Boma nett.) 
mil {3 dlog 7) 

pe dba ie ; 4 

hence ~ Beh. {5 d log BE) gly ) 


a formula effectively given by Wilson (p. 177). Let us write, following 


Mott and Jones, reel 
=| pom |’ (8) 
3eLy 
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where x stands for the expression in brackets on the right-hand side of 
eqns. (5) or (7), depending on which particular formula we use. 


§ 3. DiscUSSION OF THE CASE IN WHICH THE ELECTRONS ARE 
CoMPLETELY FREE 


In the case of completely free electrons, eqn. (7) can be used to evaluate S. 
The only unknown quantity in (7) is t which has to be known as a 
function of H. AD 

According to Wilson (p. 208), 7(k) at high temperatures, 1s given by 

7o|k|. H/C, 


where C is the ‘interaction constant’ of electrons with the lattice and is 
of the same order of magnitude as the energy of the state k. Now 
Wilson has assumed (p. 201) that C is independent of k. Therefore, in 


this case roc Bele. oe es 


Hence, according to Wilson, x=3 for all simple metals. The experimental 
values of a for alkali and noble metals are given in table 1 below, and 
x evidently varies rather widely from metal to metal. 


Table 1 
Li Na K Rb Cs Cu Ag Au 
a —6:7 2-2 3:8 2-2 0:3 =1°8 —l-] — 1-4 


On the other hand, Mott and Jones base their discussion of electron 
scattering in general on the Wigner—Seitz model, assuming specifically 
one free electron per atom and an Einstein model for the lattice. Using 
this model, let us work out the expression for 7 explicitly for any number 
of free electrons per atom—within the limitations of the Wigner-Seitz 
model—neglecting, in the first place, any distortion of the energy surfaces. 
It will be found that the variation of the number of free electrons per 
atom leads to a variation in the thermoelectric power ; and hence one 
could partly explain in this way the differing values of x. 

On this basis we have 


i ze 
==A|k| ro | f(2| k | rp sin $0)(1—cos0) sind dd . . (10) 
0 
(cf. Mott and Jones, eqn. (24), p. 254) where A is a constant independent. 


of |k| or H, ry is the atomic radius, and f(x) is given by 


(sin a—a cos a)? 
t(fjs= ee . . . . . . . (11) 
The integral above can be simplified to give 
ro |k |3/F(é), i sey you Meena eee lees 
where t= 2] k |r, 
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and LO= 5 | 


2t1]— cos « (iia, OY 


eee ceo 
0.4 — cos Z —_—-— 
Boe ea 4 


(13) 


Col or 


0 


The first integral in (13) has been fully tabulated for various values of f 
(Harvard Tables, 1949). 
From (12) it follows that 


| d log r | 2 E log zal 
d log | K | \ij=x, d log | k | \nj=n, 


a (sin t—t cos t)? 
=3— [en . . . . (14) 


where to=2korp=3:84n,1* eee 2, Ferrel) 


and nm, is the number of free electrons per atom. 

The values of (d log z/d log|k|))—;, calculated from eqn. (14) for 
various values of m , and the corresponding values of 2 are given in 
table 2 (see also fig. 1). 


Table 2 


No © | 1-2) 1-1] 1-0} 0-9} 0-8) 0-7/ 0-6) 0-5 | 0-4) 0 


(d log z/d log |k|)i=z,=a@, say | 3 | 2-6] 2-4] 1-9] 1-7] 1-5] 0-9) 0-4] 0-3 |—0-4 |—3 


a=(3+a)/2 (cf. eqn. (7)) | 3 | 2-8] 2-7] 2-5] 2-4| 2-3] 2-0]1-7]/ 1-6] 1-3] 0 


Fig. | 
x 
Ss 
fa. 
{ 
O 
O j re 3 20 


It appears that a range of values for « between ~2-8 and ~2:3 might 
be accounted for on the basis of an expected variation in the number of 
free electrons in these metals, although (cf. Mott and Jones, p. 98) on the 
extreme approximation of ‘ tight binding ’, one electron per atom would 
only correspond to as little as ~0-26 free electron/atom. However, no 
value of n which can be assumed on this model can make x greater than 
3 or less than 0. Consequently, the case of potassium for which x=3-8 
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and that of lithium for which x is negative cannot be explained on this 
basis. 


§ 4. EXTENSION TO AN ENERGY-DEPENDENCE H=const. k?-*, WHERE 
0<a<2 
In the foregoing, we have discussed completely free electrons where, 
therefore, H=(h?/2m)k?. In this section we shall consider a somewhat 
more general case where # is given by 


E=const. 2“ .; . . . 90 2 eee 


and as « varies between 0 and 2, (16) will correspond to the transition 
from absolutely free electrons (or an almost empty band) to a completely 
full band. 

Since we may always reasonably expect that the scattering due to 
lattice vibrations will only depend on the speed of the electrons,* and 
that the faster electrons will be less scattered than the slower ones 
(cf. Mott and Jones, p. 311), we may assume 


t=const.| grad H|*® .) <) .  e 


where y is a positive numerical constant depending on the number of 
free electrons per atom in the metal (cf. previous section). 
Hence from eqn. (4) we get 


o(#)=const; bY 8s te ee ce 
—a(y+1 
and hence — ee : « 2. at se 


For convenience let us now give y its limiting value of 3 in what follows. 
This corresponds exactly to Wilson’s model and sufficiently closely to 
the Wigner—Seitz model for ny =1.+ 

2(3— 2a) 


Then t= ———-. oe ee ee 


2—a 
Figure 2 shows x as a function of «. 


When the band is practically empty, i.e. when « is nearly zero, a is 
positive and takes practically its full value 3. It is obvious that 2 cannot 
be greater than -+-3, which is its limit towards the empty band side. 
On the other hand, when the band is nearly full, i.e. « is nearly equal to 2, 
w becomes negative and very large, tending towards infinity. This 


a eer eat eo ee 


* We are, of course, ignoring scattering due intrinsically to a possible 
Bragg-reflection of the electrons. Such scattering corresponds essentially to 
the Peierls’ ‘Umklapp-prozesse’ of which there is no experimental evidence in 
the resistance of the alkali and noble metals. 

+ The inclusion of a more precise value for y(0<yS3) throughout the filling 
of the band will not alter the limits for « deduced below. 
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implies that a very high positive thermoelectric power is possible when 
the band is almost full (cf. also Mott and Jones, p. 312). On the basis 
of such an argument, the case of lithium may be explained. It should, 
however, be noted that the band has to be markedly distorted by the time 
that x changes sign (H «k1/?) on this model; on the other hand, z then 
varies very rapidly with further distortion of the band. 


§ 5. INFLUENCE oF A Non-SPHERICAL FreRmMt SURFACE ON THE 
VALUE OF & 


So far we have evaluated x assuming a spherical Fermi surface. Let us 
now consider a Fermi surface which is not a sphere but can be divided 
up into segments whose surfaces are spherical or planar. These surfaces: 
may be supposed oriented in k space in such a way that cubic symmetry 
of the Fermi surface is retained, since we only wish here to consider 
metals of cubic symmetry. 


Fig. 2 


Considering then eqn. (4), it follows that o,(#)—the contribution from 

a particular ‘ segment ’ to the integral—can be written in the form 
ina Ae eg eee ey 
It can then be shown with analogous assumptions that «;<3 for every 
segment of the type considered. The total contribution will therefore 


be given by ise (M\= SAE... t. (22) 
dlogo(E)  XA,aH™ 3 
Hence “d log E = PHETS : hee SG rth Oa te ( 3) 
Since each «,<3, it is trivial to show that 
d log o(£) 
ie ee) Cn a, het, eee 
aoe G4) 


In other words, it appears that 2 ought not to exceed +3 in a metal 
of cubic symmetry, even when the Fermi surface is not itself spherical. 
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Consequently one is still unable to explain the value of for potassium, 
which is significantly greater than 3. 


§ 6. CONCLUSION 


While it appears possible on the band theory to account for the 
thermoelectric power at high temperatures of the noble and alkali metals, 
with the exception of potassium, this metal appears specifically anomalous. 
We might either suggest that in potassium the variation with energy of 
the thermal scattering is markedly different from that usually assumed 
or proposed. Alternatively, in view of earlier experiments by MacDonald 
(1952) on the electrical resistance, there may be additional band-to-band 
scattering, which would then depend on the density of states in the over- 
lapping bands and their possible variation with temperature. One 
difficulty, however, in this hypothesis is that, while Rb and Cs also 
show an anomalous resistance variation at higher temperatures, both of 
these metals appear to fit sufficiently well into the theory already discussed. 

In general, however, it should be noted that a considerable deviation 
from the idealized free-electron model is in any case called for in the 
theory to account for the range of thermoelectric power in alkali and 
noble metals. 
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ABSTRACT 


Nucleation from solid solution is investigated by means of dielectric 
measurements. The material investigated is a crystalline mixture of 
‘ equivalent ’ ketones, a system analogous to a mixture of isotopes which 
tends to segregate. During successive stages of nucleation within an 
unstable solid solution, dielectric constants and loss angles increase and 
pass through maxima. Larly stages of nucleation represent short range 
ordering and are structure insensitive, while later stages represent long 
range ordering and are comparatively structure sensitive. 


§ 1. INTRODUCTION 

THE formation and development of nuclei of new phases out of a homo- 
geneous solid solution (Smoluchowski 1951) has been studied extensively 
in alloys, but has so far not been investigated for dielectrics. However, 
certain types of polar organic materials offer exceptionally favourable 
conditions for the study of nucleation. Dielectric measurements during 
the break up of a mixture yield information not only with regard to the 
mechanism of nucleation, but also with regard to the dielectric properties 
of small aggregates. 


§ 2. THE MATERIAL UNDER INVESTIGATION 

The material investigated is a mixture of equal amounts of five 
equivalent ketones (Daniel 1950). Chemically, all these are derived from 
a saturated paraffin of given chainlength, by substitution of one oxygen 
atom for two hydrogen atoms, at a link not less than three carbon atoms 
from the end of the chain. The equivalent ketones used have a chain- 
length of eighteen carbon atoms, and can be described by the general 
formula: C,H. - CO.Ci7-m Ha;-2m. They have very closely similar 
physical and chemical properties, in consequence of their similar structure 
(Daniel 1950). The system used is a special case of a group of substituted 
long chain paraffins of equal chainlength, and other equivalent systems 
can be envisaged, with different chainlengths and different substituted 
sidegroups. 

Although equivalent ketones are very similar to one another they 
have a tendency to unmix in the solid state, because in a crystal of a 


* Based on Report L/T 291 of The Electrical Research Association. 
+ Communicated by the Authors. 
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pure ketone all keto-groups are on one level, and their electrostatic 
interaction contributes a large term to the bonding energy of the structure. 
Mixtures of equal amounts of a sufficient number of such ketones are 
thermodynamically stable either as a single-phase solid solution or a 
multi-phase mixture of the constituents, depending on the number of 
components present. The single phase form occurs, however, not only 
when it is thermodynamically stable, but also as a metastable form. 
If this occurs, it changes in due course into the multi-phase form, but the 
course of precipitation may be very prolonged. The time taken depends 
on the number of equivalent ketones in the mixture, and is the longer 
the larger the number of ketones. For the five ketones used in this 
work, completion of the phase change takes some months. 

The multi-phase form of the mixture of five ketones is thermodynamic- 
ally stable up to its melting point at 27-5°c, and consists of five separate 
phases, each made up of one ketone with a small amount of the others in 
solid solution. The crystal structure of these phases is orthorhombic, as 
for the pure ketones. 

The single-phase form of this mixture is not stable, but has the same 
melting point as the multi-phase form, within the limits of error. This 
metastable form undergoes a change of crystal structure at a temperature 
of 14:-5+1°c. It is orthorhombic, like the multi-phase form, below the 
transition temperature and hexagonal above the transition. 


§ 3. THEORETICAL CONSIDERATIONS 


Nucleation from Solid Solution 

The simplest case of nucleation in general is the formation of a solid 
or liquid phase out of vapour. This is reasonably well understood 
(Volmer 1939). The crucial point of the formation of a nucleus of a 
condensed phase is that the surface area of a solid or liquid has a higher 
free energy per molecule than its interior. In consequence small nuclei 
are comparatively less stable than large ones, and nuclei have to achieve 
a minimum size before they become thermodynamically stable. This 
has the effect that nucleation depends a great deal on circumstances 
which enable nuclei to reach the stable size, and makes nucleation, even 
in its simplest forms, erratic and difficult to control. 

Nucleation within a solid matrix, as distinct from nucleation from the 
vapour or a liquid solution, is further complicated by size relationships. 
If a nucleus has not exactly the same specific volume and shape as the 
matrix, its formation involves stresses. These affect the stability of the 
nucleus, and if nucleus formation involves a great deal of stress, it may 
become practically impossible in the interior of the matrix. In this case, 
nucleation depends very much on the texture of the matrix so Gade 
‘structure sensitive ’ because nuclei can usually arise much more easily 
on surfaces or other disturbances of the matrix than in the interior. 

Theoretically, structure sensitive nucleation is exceedingly complex, 
and the existing theoretical work (Smoluchowski 1951, 1952) is restricted 
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to ‘structure insensitive ’ nucleation and the erection of nuclei within 
the solid matrix. Most of the existing theories use a particular model 
of a solid solution for which certain simplifying assumptions are made. 
This model is a binary mixture of two substances, A and B, of equal 
crystal structure, for which the tendency to form two phases can be 
expressed by a single quantity : 


V= Viap—3(Vaat Vo) . ° . . . . . (1) 


where V,,, V,, and Vy, are respective energies of pairs of atoms (or 
molecules). With this assumption the free energy per gramme atom of 
mixture of concentration « can be written : 

S(«)=VZa(l1—a)+kT'(« log «+(1—a«) log (1—«)) aul ieeeec (2) 
where Z is the number of nearest neighbours in the structure. Equation (2) 
implies that for negative V the material is single-phase at all temperatures. 
while for positive V it breaks up into two phases below a temperature 
which depends on the concentration. 

In discussing nucleation from solid solution Smoluchowski (1951) takes 
as his starting point the concept of random fluctuations of concentration 
which exist in any mixture. He shows that where the mixture is single- 
phase in equilibrium, positive V, that is a tendency towards precipitation, 
increases the amplitude of random fluctuations. In other words, the 
single-phase material contains very loose and indefinite aggregates of the 
two constituents, which are the more distinct the closer the material is 
to two phase conditions, but which are continuously formed and 
redissolved. If the material is two phase in equilibrium, aggregates of 
this kind act as embryonic nuclei which may grow into stable nuclei if 
the free energy conditions permit it. 

The experimental evidence on nucleation consists mainly of investiga- 
tions of alloy systems (Geisler 1951). It confirms the theoretical ideas 
qualitatively, but is rarely suitable for quantitative tests of the theories, 
because each alloy system tends to introduce its own complications 
(Smoluchowski 1952). However, x-ray evidence of the earliest stages 
of nucleation (Guinier 1952) that are accessible shows that early nuclei 
have structures which are neither quite like that of the final precipitated 
phase, nor quite like that of the matrix. This confirms the ideas of 
embryonic nuclei arising from random fluctuations within the single-phase 
matrix. 

The mixture of five equivalent ketones is closely related to the idealized 
model of a binary mixture with a solid solubility gap, the properties of 
which are represented in eqn. (2). In fact, a binary system of two 
equivalent ketones is closer to this model than any alloy system can be, 
because of the great similarity of equivalent ketones. 

A mixture of five equivalent ketones is not much more complicated. 
than a mixture of two ketones, because its free energy conditions can 
still be represented with the help of a single constant, 


Ves Vip—3(Vaat Von), 
4U2 
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where A represents a particular equivalent ketone and B any other 
equivalent ketone. The quantity V can in principle be deduced from 
measured properties of equivalent ketones (Daniel 1950), but the available 
data are at present not quite sufficient. ’ 

The characteristics of nucleation in the mixture of five ketones are 
clearly a matter for experiment, but certain points can be made a priori, 
since the ketones are equivalent and present in equal proportion. Each 
nucleus must consist predominantly of one ketone, with a proportion of 
the other four in solid solution. If any ketone forms a nucleus, it thereby — 
supersaturates the solid solution with regard to the others so that groups 
arise of five nuclei. Such groups may be separated by regions of undis- 
turbed solid solution. However, once a stable group has been formed, 
the growth of its nuclei propagates the nucleation, since each growing 
nucleus supersaturates the solution out of which it grows. 


Dielectric Properties during Nucleation 
In the present investigation the process of nucleation is followed 
largely by means of dielectric measurements, and it will be necessary to 
consider the significance of the dielectric data. 
In general, the static dielectric constant of a material can be expressed 
by an equation (Frohlich 1949) : 


3e, 47N mm* 
2etn? 3 kT 


2 
€-— n= 


(3) 


where n? is the square of the refractive index and N the number of 
polarisable units (e.g. dipoles) per unit volume. The expression mm* 
has the dimensions of a dipole moment, squared, and depends on the 
configuration of dipoles as well as on their moment when in isolation. 
Equation (3) is derived under the assumption that the structure of the 
dielectric is uniform within a spherical region small compared with the 
dimensions of the test capacitor, i.e. in the sense of this expression a 
material is homogeneous if its structure has inhomogeneities only within 
dimensions that are small, as compared to 10-3 or 10-4 em (the thickness 
of dielectric used being 5x 10-2 em). This means that nucleating material 
can be considered as homogeneous in the sense of eqn. (3) if there are 
many nuclei within dimensions of the order 10~4-10-%, that is, if average 
nuclei are smaller than about 10-°-10-5 cm. In that case the material 
is indeed no less homogeneous than a protein with large molecules. If 
aggregates are larger, the dielectric properties of the material as a whole 
would have to be treated according to the equations applying to mixtures 
of dielectrics with different permittivity (Rushman and Strivens 1947). 

If nucleation is structure sensitive, and groups of nuclei occupy an 
unknown and more or less small fraction of the material. dielectric 
measurements can yield little significant information about the condition 
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of nuclei. All that can be said in this case is that the results of measure- 
ments will vary greatly with the texture of the single-phase solid solution, 
so that they are likely to be badly reproducible. 

If nucleation is structure insensitive, the material in the course of 
nucleation consists of a mass of nuclei with boundary regions between 
them. The size of nuclei may obviously vary in space and time, but 
they may be expected to contain initially no more than some tens to 
thousands of molecules each. This follows from the metallurgical 
evidence on early nuclei (Guinier 1952) as well as from first principles, 
since diffusion over large distances is unlikely in the early stages of 
nucleation. ‘This means that early nuclei are likely to be smaller than 
10° ®-10~° em, so that the material as a whole is homogeneous in the sense 
of eqn. (3). In this case, the dielectric constant measures mm* for an 
average nucleus, where a nucleus is taken to include an appropriate part 
of the boundary region dividing it from the neighbouring nuclei. For 
larger nuclei it may be necessary to consider nuclei and boundaries 
separately as regions of different dielectric constant (Rushman and 
Strivens 1947). 

The quantity mm*, which can be deduced from measurements on a 
homogeneous material, is a measure of the tendency of dipoles to be 
parallel to each other (Fréhlich 1949). This quantity is large if dipoles 
tend to be parallel, small if they tend to be antiparallel. As a rule the 
tendency towards parallelism or otherwise is determined by the force 
fields of neighbouring molecules, i.e. short range forces. However, in 
the case of ketones dipolar interaction, which has a long range, plays an 
important role in bonding, and so both long and short range forces come 
into play. 

The loss angle gives information, as it were, from the viewpoint of the 
individual dipole. According to Debye, the loss angle is given by an 
equation : 


(€,—n?)wr 


Oh eareer er 
eo ecoeT 


(4) 
while the simultaneously measured real dielectric constant ¢, is given 
by 

€,—n 


€—V= 
1+ w?7? 


where n? is the square of the refractive index, w/27 is the frequency in 
cycles/second and 7 the relaxation time. 7 is given by: 


moaexph— HiT sg hs. (8) 


where A is approximately constant, and # represents an energy barrier 
between two positions of a dipole. For the purposes of eqn. (6) it is 
assumed that H is the same for all dipoles in the material. This in 
in general not true, in that a large number of values of H may exist if 
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all dipoles are not surrounded by identical force fields. In the case of 
precipitating aggregates we may expect a distribution of values of H 
and thus a distribution of relaxation times, with correspondingly modified 
eqns. (4) and (5). The values of # in the distribution of relaxation times 
give information about the various force fields surrounding individual 
dipoles. 


§ 4. EXPERIMENTAL PROCEDURE 


Nucleation in the mixture of five ketones was investigated using two 
experimental methods: dielectric measurements, and X-ray powder 
photographs. 

Dielectric measurements were carried out using the cooling capacitor 
and metal desiccator described elsewhere (Daniel and Stark 1952). With 
this apparatus the specimen is kept free from contamination and may be 
heated and cooled during the measurements. The arrangement also 
permits quenching by passing cold fluid through the hollow bottom 
electrode of the capacitor. An accurate Schering bridge was used 
(Daniel and Stark 1952), at frequencies of 50 c/s, 800 c/s and 50 ke/s. 

The process under investigation is a change of dielectric properties 
with time and temperature, and it was found impracticable to follow it 
simultaneously at all three frequencies. For this reason, heat treatments 
had to be repeated three times to obtain results at all frequencies used. 
This procedure is justified by the reproducibility of the results. However, 
it makes important a certain source of error which has been already 
mentioned (Daniel and Stark 1952). When a long-chain material is 
solidified in a temperature gradient it becomes preferentially orientated. 

Preferential orientation markedly affects the dipolar part of the ~ 
permittivity of a ketone, since dipoles reverse by rotation or twisting of 
the long-chain molecules. With the present material quenching produces 
very marked preferential orientation, and even slow cooling of the 
desiccator box does not quite eliminate it. Fortunately, the dielectric 
constant of the material in its hexagonal form can be used as a measure 
of preferential orientation : the higher its value, the more the material 
is oriented. If the dielectric constant of the hexagonal form is measured 
prior to a heat treatment leading to nucleation, the effects of preferential 
orientation can be eliminated, or allowed for by a numerical correction. 
Since this correction is mostly not large and is somewhat tedious, it is 
applied only where necessary for purposes of the argument. In general, 
uncorrected results are given. The details of the corrections applied 
are described in E.R.A. report L/T 291, but are omitted here for the sake 
of brevity. 
fe eh earn ee a a eke was carried out using an apparatus 
arn ¥ aa ze Be 2 allows photographs to be taken at tempera- 

emperature, and where six successive powder 
photographs can be recorded on the same film. 
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§ 5. EXPERIMENTAL RESULTS 
Dielectric Evidence 


The type of result obtained by the dielectric measurements is shown in 
fig. 1, which shows the dielectric constants and loss angles obtained when 
the material is cooled and subsequently heated again, after varying 
periods at a low temperature. The measurements recorded were carried 
out at 800 c/s. The figure illustrates clearly the complicated effects of 
heat treatment. On cooling from the melt, the dielectric constant is 
first constant, then rises to a maximum, and finally decreases, when the 
specimen is kept at a low temperature. On subsequent heating, it rises 
again and passes through a maximum, which is the less pronounced, and 
occurs at the higher temperature, the longer the specimen had been kept 
at a low temperature. 
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The dielectric constant of mixed ketones on cooling, and heating atter ageing 
below 4°c. 


The phenomena observed depend in a complicated way on time and 
temperature. In order to study them in more detail, heat treatments 
had to be devised which could be repeated conveniently. The most 
useful of these was obtained when the specimen was cooled very slowly, 
in a refrigerator set finally to about 6°c. Figure 2 shows the results of 
this heat treatment, at all three frequencies used, for both dielectric 
constant and loss angle. The results are given as a function of the time 
of cooling, the temperatures reached at a given time being marked on 
the time scale. 

The phenomena observed are best considered in several stages. How- 
ever, before doing this it will be necessary to record the x-ray evidence. 
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Fig. 2 
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Crystallographic Evidence 


The heat treatments which give large dielectric effects as recorded in 
figs. 1 and 2, have only very small effects on the crystal structure, even 
though there are appreciable differences between the structures of the 
single-phase solid solution and the stable multi-phase form. 

For the mixture of five ketones changes of crystal structure consist 
in an expansion or distortion of the (ab) plane, at right angles to the 
¢ axis (Daniel 1950). The area ab of the basal plane (ab) of the unit 
cell can be taken as a measure of the specific volume, while the shape of 
the unit cell may be characterized by the angle 40, which measures the 
deviation of the basal plane from hexagonal shape. For a hexagonal 
structure 06=60°. 


Figure 3 shows the cross sectional area ab and the angle 40 for the’ 
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mixture in its various forms and for a pure equivalent ketone. The 
single-phase form is hexagonal above 14-5--1°c and orthorhombic below 
this temperature (as shown by the angle $0). The area:a x 6 is continuous 
through this transition, which seems curious. However, a small volume 
change at the transition cannot be excluded on the basis of the available 
evidence. Between about —5°c and the melting point the single-phase 
form has a very high thermal expansion, as compared with that of the 
pure ketone. he 
The mixture approaches stability only very gradually, so that it is 
difficult to decide when stability has been reached. However, the 
specimen whose volume is recorded in fig. 4 was aged for five months at a 
temperature of approximately 5°c, and may be taken to be near to 
stability. The stable mixture shown has the same thermal expansion 
as the pure ketone, but a somewhat larger volume. This implies that its 
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bonding is similar to that of the pure ketone, but that the material 
contains some of the other ketones in solid solution as is to be expected. 

The early stages of nucleation, though exhibiting large dielectric effects, 
give x-ray powder photographs very similar to those of the single-phase 
form. Figure 4 shows the changes of structure corresponding to a heat 
treatment very similar to that which gave the dielectric effects recorded 
in fig. 2. For both heat treatments recorded the specimen had been 
cooled more or less gradually to 7-5°c, kept at that temperature until 
3-31 hours and then heated again. The ageing at 7-5°c produces a slight 
contraction which is undone on reheating. 


Fig. 4 
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Structure Insensitive Nucleation 

The X-ray record may be used to decide whether or not the changes 
on ageing occur uniformly throughout the material, i.e. whether the 
observed nucleation is structure insensitive. To this end, the X-ray 
lines 110 and 200 were photometred, to ascertain whether they broaden 
as they shift; a Hilger microphotometer was used. Details of the 
microphotometer record, as well as a reproduction of a typical set of 
photographs, are shown in E.R.A. report L/T 291. The maximum line 
shift observed (for 200) is only 0-05 cm, about half the width of the line, 


Nucleation in a Dielectric 1381 


so that the data cannot be evaluated very accurately. However, the 
fact that no change of line width is observed while the lines shift may be 
taken to mean that changes occur simultaneously in at least two thirds 
of the bulk of each of the specimens examined. This implies that the 
changes observed are not restricted to surfaces or other anomalous parts 
of the specimen. Nucleation is essentially structure insensitive. 

The crystallographic evidence is confirmed by a result of the dielectric 
measurements which will be discussed in more detail later. The highest 
value of the dielectric constant obtained on cooling during heat treatments 
as shown in figs. 1 and 2 was found to be reproducible to better than 
+6%, once allowance is made for preferential orientation. This result, 
confirmed for eleven separate heat treatments, would be very unlikely 
if nucleation were structure sensitive. The small scatter of the results 
suggests that nucleation occurs simultaneously in the larger part of the 
specimen, with a few exceptions in some of the crystals of the matrix. 


The Successive Stages of Nucleation 


The fact that nucleation proceeds uniformly within the bulk of the 
material makes it possible to interpret the observed changes of dielectric 
properties. These changes correspond to a process in the material, and 
the dielectric properties at each stage represent average properties of nuclei. 

The process of nucleation appears to be divisible into a sequence of 
more or less separate stages which succeed one another as the temperature 
gradually decreases. These stages will be discussed in turn. 


The Hexagonal Form 


When the specimen is single-phase in the hexagonal form, no nucleation 
occurs within times of weeks or months. This is not surprising, since 
the multi-phase form has an orthorhombic crystal structure and nucleation 
within the hexagonal form would presumably produce large stresses. 

The dielectric constant of the hexagonal form is almost independent 
of temperature. It is lower than that of the liquid mixture, but higher 

_than that of the pure ketone, for which the dipolar contribution is zero. 
This means that the dielectric constant of the hexagonal form contains 
a comparatively small dipolar contribution. It is curious that this 
contribution is not larger and that it does not cause the dielectric constant 
to rise with falling temperature, as it should according to eqn. (3), if 
mm* is constant. No convincing explanation is available of this fact, 
but it seems that in the hexagonal form dipoles of neighbouring different 
ketones tend to oppose each other, so giving a comparatively small mm* 
which decreases with temperature. 


The Earliest Stage of Nucleation 


On cooling the single-phase specimen through the transition from 
hexagonal to orthorhombic, the dielectric constant rises suddenly. The 
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sudden rise at 13°c in figs. 1 and 2 is reproducible, and corresponds to 
the sudden onset of the structural transition when the material is cooled 
to the lower limit of its temperature hysteresis. When the change is 
complete the dielectric constant rises slowly with falling temperature, 
while the loss angle at all three frequencies remains as low as for the 
hexagonal form until the temperature falls to about 10°c. On further 
cooling the loss angle also begins to rise. 

Between the temperature of the transition and about 10—-11°c the 
dielectric constant is almost independent of time. This is demonstrated 
by a temperature cycle (fig. 5) when the specimen was cooled to 11°c 
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and subsequently heated, the cycle taking 5 hours to complete. A small - 
permanent change occurs at the lower end of the temperature range, and 
this is connected with a very slight contraction of the material. However 
in first approximation the dielectric constant and volume are foriotians 
of temperature, but not of time. 

The time independent rise of the dielectric constant just below the 
transition temperature may in principle be explained in two ways: it 
may be a property of the unchanged single-phase form or it may be due 
to a rearrangement of molecules which depends on temperature but is so 
rapid as to appear time independent. Two arguments speak for the 
latter alternative. Firstly, the slope of the «(7’) curve indicates an 
increase of mm* with falling temperature ; for constant mm* we should 
expect a slope de/d7'=5 x 10-3, while the measured value is always larger 
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than de/dT=4X10-*. This indicates a tendency of the dipoles to become 
parallel to one another with decrease in temperature, as is to be expected 
if clusters of the individual ketones are formed. Secondly, the phenomena 
occurring between the transition temperature and 11°c appear to be a 
necessary preliminary of the more striking nucleation processes occurring 
at lower temperatures. If a specimen is quenched through this range in 
a matter of seconds, the (corrected) maximum value of the dielectric 
constant is low, e=4-1 at 800c/s on the average of two occasions, as 
against «-=4-6+0-2 when the specimen is cooled more slowly, spending 
some 5-10 minutes between 13°c and 11°c. 

The observed features of the first stage of nucleation suggest that 
it consists in an increase of short range order, in that the individual 
ketones form embryonic nuclei. The distinctness or number of these 
nuclei increases with falling temperature, and it appears that the character 
of nuclei is in a metastable equilibrium with temperature. This observa- 
tion is in agreement with recent theoretical work by Dehlinger (1951), 
who considers structure insensitive nucleation in a material where the 
formation of aggregates creates stresses. He postulates that the stresses 
increase with increasing aggregate size, thus creating a limiting size 
beyond which aggregates cannot grow. The average size of nuclei should 
be in a metastable equilibrium with temperature. 

In the mixed ketones the metastable equilibrium is restricted to a 
temperature range of only some 2°c. This is in good agreement with 
Dehlinger’s ideas, because in the temperature region in question the 
volume difference between the stable multi-phase and the single-phase 
form decreases rapidly with falling temperature (fig. 3). Stresses may 
be expected to be the larger the larger this volume difference, so that the 
limiting size should increase rapidly as the temperature decreases. 

The metastable equilibrium observed implies that the embryonic nuclei 
in question are formed and destroyed rapidly. This has the consequence 
that after a short time—say a few minutes—at a given temperature the 
material should be saturated with embryonic nuclei of the number and 
size appropriate for the temperature. 


The Stage of Maximum Dielectric Constant 


On cooling below about 11°c the character of nucleation changes, as 
shown in fig. 2. The dielectric constant rises more rapidly, and the loss 
angle begins to increase also. The increase in loss angle occurs the sooner 
the higher the frequency, and this, according to eqn. (4), indicates a 
gradual increase of the average relaxation time. The process observed 
is no longer independent of time, as shown by fig. 6, which gives the 
dielectric properties at 800 c/s for a rapid and a slow cooling sequence, 
For both heat treatments the dielectric constant reaches a maximum 
followed by a decline, but for rapid cooling the maximum is reached 
sooner, and at a lower temperature, than for slow cooling. 
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The experimental data indicate that the embryonic nuclei formed 
above 11°c become able to develop at lower temperatures and that their 
development needs time. The behaviour of the loss angle as a function 
of frequency indicates that this development consists, at least to some 
extent, in a growth of nuclei, in that molecules within them become more 
firmly bonded. 

The development of nuclei always proceeds in the same way even if 
its speed varies with temperature. This is shown by the fact that the 
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maximum dielectric constant observed during a heat treatment is, for a 
given frequency, independent of the speed of the treatment, once allow- 
ance is made for preferential orientation. Table 1 gives the characteristics 
of the condition of peak dielectric constant for eleven separate heat 
treatments. Preferential orientation is eliminated by calculating mm* 
peak/mm* hexagonal. The loss angle, which is also corrected for prefer- 
ential orientation, refers to the condition of peak dielectric constant. 
The characteristics of the heat treatments are given by the temperature 
at which the peak occurred and the time between the beginning of nuclea- 
tion at about 13°c and the occurrence of the peak. 

Table | shows that mm* ,/mm*y is constant independent of the length 
of the heat treatment, the temperature at which the peak occurs, and of 
preferential orientation. This indicates that the condition of maximum 
dielectric constant represents a certain texture through which the material 


passes in the course of nucleation. It appears to correspond to a certain 
stage in the development of the nuclei. 
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Table 1. Particulars of the Condition of Peak Dielectric Constant, 
Measured at 800 c/s 


—— 2 
' mm*,, tan 8x10 is Time 

€pt «at |—— (°0) oa Remarks 
mm” «| Uncorrected| Corrected ») | (minutes) 

5:02 | 3-54 | 1-91 1-2 1-03 1-4 30 

5:81 | 3-84 | 1-98 1-4 1-02 2-0 50 

4:56 | 3:33 | 1-86 0-8 0:77 2-8 50 

7 14.9% | 1-80 1-4 0-81 3-5 69 {| Quenched | 

from melt 

DIS DO) 2°02 1:3 1-05 ee 70 

4-44 | 3-34 | 1-79 0-7 0:72 6-1 7) 

5:88 | 3-82 | 2-02 1:7 1-22 2-3 80 

4:94 | 3-50 | 1-89 1-0 0-88 4-5 100 

5:04 | 3°52 | 1-92 1-2 1:03 4-8 110 

4-81 | 3-52 | 1-79 1-0 0-91 7-2 180 See fig. 2 

Aeay |) Sst ele) 0-8 0-83 4-8 300 
1-90 0-98 x 10-# Averages 


+ «,—the maximum dielectric constant reached during the heat treatment. 
{ e,—the dielectric constant of the hexagonal form, measured prior to the 
heat treatment. 


Table 2. Particulars of the Condition of Peak Dielectric Constant, 
Measured at 50 c/s and 50 ke/s 


mm* 
= tan 8 (corrected) | Z'y (°c) Time (minutes) Remarks 
mm* 4 
50 c/s: 
1-81 4:1x10-% 7-2 180 See fig. 2 
1-94 alex 104 4-8 75 
1-88 3°6 x 10-* Averages 
50 ke/s 
1-65 2:4 10-2 75 160 See fig. 2 
1-69 1-9) 10-* 3:7 55 
1-67 2:2 10-2 Averages 
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Although the value of mm* ,/mm*y does not significantly depend Os 
any of the recorded variables, there still exists a scatter of about +6% 
about the average, indicating an influence of chance on the progress of 
nucleation. This is to be expected for a nucleation phenomenon, even 
if it is essentially structure insensitive. , 

The loss angle associated with the peak dielectric constant is also 
constant when corrected for preferential orientation. The scatter of 
values is higher than for the dielectric constant, being +20%. However, 
the accuracy of measurement is lower in this case, because the loss angle 
usually rises rapidly while the dielectric constant is at its peak, 

The peak values of the permittivity, and the corresponding loss angles 
were measured algo at 50 c/s and 50 ke/s, giving results similar to those 
at 800 c/s, as shown in table 2. 

Considerations of fig. 3 and of the two tables shows that the maximum 
permittivity at 50 c/s and 800 ¢/s corresponds to a distinct texture, and 
the crystallographic evidence (fig. 5) shows that this condition is not far 
removed from the single phase state. For 50 ke/s, the maximum is 
lower and reached earlier, indicating that some dipoles contributing to 
the dielectric constants at the lower frequencies are unable to do so at 
the higher one ; this is understandable on the basis of eqn. (5) and implies 
a ‘ freezing in’ of dipoles with comparatively long relaxation times. 

The occurrence of a maximum of the dielectric constant at some stage 
‘of nucleation is not difficult to understand. In a pure ketone, and there- 
fore also in a stable nucleus of the multi-phase form, dipoles are arranged 
in ‘dipolar planes’, within which they are parallel to one another. 
Frohlich (1946) shows that this configuration may, in certain critical 
conditions, lead to high values of mm* and the dielectric constant. 
However, experiment (Daniel and Stark 1951) shows that pure ketones 
have a negligible dipolar contribution to the dielectric constant, pre- 
sumably because the critical conditions are never reached, dipoles being 
too firmly bound. The dielectric constant of the stable multi-phase 
form, consisting of more or less pure ketones is similarly low (fig. 1). 
On the other hand, the single-phase form also has a permittivity low 
compared with that of the liquid, due presumably to a tendency of 
neighbouring dipoles of different ketones to be antiparallel. It seems 
plausible that at some stage between the extremes of mobile but more or 
less antiparallel, and largely parallel but immobile dipoles there exists 
a condition when dipoles are mobile and tend to be parallel, thus leading 
to a comparatively high value of the dielectric constant. 

While the existence of a maximum of the dielectric constant at some 
stage of nucleation can be explained without difficulty, its observed 
magnitude can so far not be interpreted theoretically ; the comparatively 
low value observed may be due to a number of causes. Further progress 
would demand a knowledge of the structure of nuclei and the disordered 
regions between them, as well as a more detailed knowledge of the distri- 
bution of relaxation times. It is possible that dipoles in the interior of 
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nuclei have very long relaxation times, and do not contribute to the 
dielectric constant even at low frequencies. 

Information on the structure of nuclei might be obtained with the aid 
of x-ray methods analogous to those employed for alloys (Geisler 1951, 
Guinier 1952), but such work would demand single crystals and its 
chances of success are difficult to assess. However, the determination 
of a distribution of relaxation times as a function of aggregate size would 
be of considerable general interest with regard to the theory of viscous 
liquids and of plastics, which often show wide distribution of relaxation 
times and are often supposed to contain semi-crystalline aggregates. 

The mixed ketones investigated in the present report may not constitute 
the most interesting system of equivalent substances for the investigation 
of dielectric properties of small aggregates. Long-chain secondary 
alcohols might be more interesting in this respect, for two reasons : 
because they exhibit comparatively high dielectric constants (Meakins 
and Welsh 1951), and because a theory (Sack 1952) exists for the dielectric 
constant and loss angle of chains of linked hydrogen bonds, such as would 
be formed in small aggregates of secondary alcohols. If such alcohols 
show nucleation similar to that in mixed ketones they may prove a 
fruitful system for investigation. 


Later Stages of Nucleation 


The stage of maximum permittivity is well at the beginning of the 
process of nucleation, and is followed by a long sequence of later stages. 
The dielectric constant and loss angle fall monotonously as these stages 
are traversed at constant temperature, and give little information on 
their nature. However, information can be obtained by reheating the 
specimen. As fig. 1 shows, heating after ageing at low temperatures 
reverses the development of nucleation, and the specimen appears to 
traverse again the condition of maximum dielectric constant. The 
conditions in which the peak occurs on reheating give information on the 
stage of precipitation that had been reached. 

Figure 7 shows that if a specimen has been aged only for a short time 
after passing the peak of the permittivity, it develops a peak on reheating 
which is somewhat lower than the peak on cooling. This peak occurs at 
a temperature which is not very much higher than the temperature of 
the original peak. Figure 1 shows that the peak on reheating occurs at 
higher temperatures and becomes reduced in magnitude as longer ageing 
precedes reheating. 

X-ray evidence on reheating after short ageing, shown in fig. 4, confirms 
the results of dielectric measurements in that the changes on cooling are 
undone by heating. For longer times, the changes are reversed at a 
higher temperature, again in accordance with dielectric evidence. ; The 
x-ray data also show that the specimen is no longer uniform when it has 
been aged for a matter of days or weeks: some parts of it retain their 
low temperature form up to a higher temperature than do others. 
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The experimental evidence on later stages of nucleation can be inter- 
preted from the viewpoint of thermodynamics and from the viewpoint 
of kinetics. 

The thermodynamics of the different stages of nucleation can be 
represented very simply for a mixture of equivalent ketones, because 
all constituents have identical thermodynamic properties. In any stage 
of nucleation, which occurs uniformly in the material, it is possible to 
define an average free energy G per mole of one (any one) of the ketones, 
and this quantity describes the stability of the material as a whole. 


Fig. 7 
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60 minutes 
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Dielectric constant € at 800c/s 


Loss angle tan § at 800c/s (dash line curve) 


Temperature , °C 


Cooling followed by heating. 


Figure 8 shows a schematic diagram of the free energy conditions as 
revealed by the experiment. Graph 1 represents the liquid, graph 2 
the single-phase form, graph 5 the stable multi-phase form, while graphs 
3 and 4 are characteristic of two successive stages of nucleation. The 
figure implies that the intermediate stages of nucleation are more stable 
than the single-phase form at low temperatures, but that above some 
a aa 7’, and 7’, the intermediate stages revert to the single-phase 

The free energy equals G—=H—T'S, where H is the total energy and 
S the entropy. Hence the slope of the graphs is (roughly) indicative of 
S and their height of H. The explanation of the observed data demands 
a comparatively large increase of the total energy of the intermediate 
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forms as compared with the stable form. This may be attributed to the 
energy of boundary regions. 

As for the kinetics of the later stages of precipitation, the most salient 
fact is that the development of nuclei becomes less uniform the further 
it proceeds. The later stages of precipitation are increasingly structure 
sensitive. This is probably due to the breaking up of the coherence of 
the single-phase matrix by the formation of compact nuclei of the new 
phases. Boundaries between them are probably no longer regions of 
gradual transition, but clear cut interfaces across which diffusion can no 
longer proceed as simply as before, so that further growth depends on 
the configuration of the interfaces and so on chance. The formation 
of new interfaces is suggested by the appearance of the mixture in the 
later stages : while the material is transparent up to the stage of maximum 


permittivity, or a little beyond, it becomes opaque as nucleation proceeds 
further. 


Fig. 8 
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Schematic representation of the free energy conditions. 


Once stable nuclei are formed, their further development is a process of 
grain growth and of uphill diffusion, i.e. diffusion against the concentration 
gradient. Uphill diffusion in the mixed ketones might perhaps be studied 
by methods analogous to those employed in a two-phase alloy (Daniel 
1948), namely by comparing the speed with which corresponding stages 


Axe 
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are reached at different temperatures. In favourable circumstances such 
work would allow the measurement of coefficients of diffusion, or at least 
of their activation energy. This would be of particular interest in mixed 
equivalent ketones which may be assumed to have equal coefficients of 
diffusion, so that one would measure effectively coefficients of self diffusion 
for different conditions of texture. 


§ 6. CONCLUSIONS 


As regards nucleation, the chief result of the present work is that the 
early stages of nucleation in the mixed ketones represent short range 
ordering and are structure insensitive, while the later stages represent 
long range ordering and are comparatively structure sensitive. 

The earliest stages of nucleation are of particular interest. In the 
temperature range above about 11°c the material appears to contain 
large fluctuations of the concentration which are in a condition of meta- 
stable equilibrium : their magnitude (or number) depends on temperature, 
but not significantly on time. Cooling below 11°c brings about a change 
in the character of nucleation, in that nuclei are now able to grow in time. 
This gradual, but fairly sharp change seems to mark a transition from 
embryonic to stable nuclei. 

Once stable nuclei are formed, they soon make up the bulk of the 
material and later stages consist of a gradual development of the nuclei 
from loose aggregates to crystals of the fairly pure ketones. This 
development is reproducible in its early stages, in that the material 
always traverses the same sequence of textures. 

The dielectric properties of the different stages of texture can so far 
not be interpreted in detail. Further experiments are needed to deter- 
mine the size and structure of nuclei, and the distribution of relaxation 
times. The results of such work may be of general interest for the theory 
of viscous liquids and plastics, where long relaxation times may be due 
to the presence of semi-crystalline aggregates. 

The methods used in the present investigation need not be restricted 
to the mixture of five ketones used, since this system is only one of several 
systems of equivalent substituted paraffins that can be envisaged. Of 
these, equivalent secondary alcohols would be of particular interest, in 
view of results which indicate high dielectric constants. 

Systems of equivalent substances constitute rather cumbersome 
research tools because of the chemical work involved in their preparation. 
However, they offer the advantage that conditions are thermodynamically 
simple, and that they can be altered in small steps, by altering the chain- 
length, the nature of the side groups, or the number of equivalent 
substances. The use of such systems may be worthwhile for the 
systematic investigation of nucleation, diffusion, and of the physical 
properties of small aggregates. 
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SUMMARY 


The unit cell of an «-phase in the Al-Fe-Si system is found to be 
hexagonal with a=12:3 A and cy=26-2 A, and differs from the cells 
previously reported for other phases of closely similar composition. 
The most prominent Brillouin Zone associated with the structure of this 
compound is compared with that of one of these other phases with special 
reference to the possible ‘ effective valency ’ of transitional metal atoms 
in electron-rich environments. 


§ 1. INTRODUCTION 


THIS paper presents the results of an x-ray examination of single crystals 
of the «-phase in the Al-Fe-Si system prepared by Pratt and Raynor 
(1951 a). A number of other workers have also examined the aluminium- 
rich region of this system. All are agreed on the occurrence of at least 
two ternary phases, usually denoted «(AlFeSi) and f(AlFeSi), although 
Phragmén (1950) has called them c- and m- respectively, and Armand 
(1952) has suggested that the «-region contains three distinct phases 
denoted «,, «) and «3 in order of increasing silicon content ; «; may be 
outside the range of alloys examined by other workers. 


The results of analyses of extracted crystals of these compounds are 
summarized below : 


Fe (wgt.%) Si (wgt.%) 
Dix and Heath (1928) o 30 8 
Phragmén (1950) c 31-9 5:6 
Pratt and Raynor (1951 a) « 32-1-32-7 8-4-10°3 
Armand (1952) Oy 27:3 7:0 
XL 29-2 11:3 
Xs 35:3 12:8 


The unit cell of c-(AlFeSi) has been reported by Phragmén (1950) to 
be cubic with ayj=12-523 kx, and the structure is considered to be isomor- 
phous with that of c-(AlMnSi), for which ay= 12-625 kx. Phillips and 
Varley (1943) have reported that in the quaternary Al-Fe—Mn-Si system 
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these two phases are the end members of a continuous series of solid 
solutions. Armand’s x-ray examination (1952) confirmed that his O4- 
phase corresponds to the cubic c-(AlFeSi). His «-phase, on the other 
hand, has a hexagonal unit cell whose a- and c-parameters are multiples 
of 4:95 kx and 7-07 kx respectively. 


§ 2. RESULTS OF THE x-Ray EXAMINATION 


We have been able to make an x-ray examination of single crystals of 
all the samples of «(AlFeSi) described by Pratt and Raynor (1951 a) 
and find none of them to be cubic. They are, in fact, hexagonal, and 
the cell dimensions, as measured by layer line measurements on oscillation 
photographs, are 

@yg—12:3+0-1 A, Cy=26-2+-0-2 A. 

The Laue symmetry and systematic absences are consistent with space 
group D§,—P6;/mmc, although no test has been made for the presence 
or absence of a centre ofsymmetry. Additional evidence for the hexagonal 
symmetry of «(AlFeSi) may be seen in the external form of the crystals 
shown in fig. 18 in the paper by Pratt and Raynor (195la). There 
appears to be no simple relationship between the hexagonal unit cell 
found for Pratt and Raynor’s crystals and that given by Armand for his 
%»-phase. 

Rotation and oscillation photographs of single crystals of «(AlFeSi) 
of composition 32-5 wgt.% Fe and 8-7 wgt.% Si, taken with Mo Ka 
radiation, show a marked grouping of very intense reflections at Bragg 
angles, 20, close to 20° and no other intense reflections with 20 less than 33°. 
The planes (Akil), listed below, which give rise to these reflections may 
therefore be considered to define a very prominent Brillouin Zone in wave 
number (k) space, with faces parallel to the planes (Aki) at a distance 
p=1/2d(hkil) A~! from the origin of the zone, d(hkil) being the interplanar 
spacing. 


(hkil) oe ae p(recip. A) 
(000-12) 2 0-229 
(202-11) 12 0-230 
(303-10) 12 0-238 
(3257) 24 On2ad 
(3360) 6 Cree 


This is the only prominent zone with faces closer than 0:40 A-! to the 
origin ; it is shown in fig. 1. 


§ 3. DIscuUSSION 


It is clear that there are at least two ternary compounds in the a-region 
of the Al-Fe-Si system, that having the lower silicon content being 
cubic (c- of Phragmén, «1 of Armand) and that having the higher silicon 
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content having a hexagonal structure (« of Pratt and Raynor). i 
the latter is really different from the hexagonal structure of Armand’s 
a5-phase, and if his «;-phase is also included in this region, there must be 
four different ternary compounds in the ‘« ’ region of the system, and 
further careful metallurgical work will be required to determine their 
precise relationship. 


The first prominent Brillouin Zone of «(AlFeSi). 


It would be interesting to consider the influence of their structures 
on the stability of these compounds, but sufficient data are not available. 
By considering the composition ranges of «(AlMnSi) and «(AlFeSi) 
Pratt and Raynor (1951 b) have suggested that both phases are stabilized 
by an electron/atom ratio of about two electrons per atom. If, as seems 
very likely, these compounds may be regarded as electron compounds, 
we can at least compare the stability of two of them (c- and «) on the 
basis of the x-ray measurements already made. 

For «(AlFeSi), the very prominent Brillouin Zone described in the 
previous section has a volume (V) 0:06444-%. We might expect this 
hexagonal phase to be stabilized at a composition having a ‘ free’ 
electron/atom ratio slightly less than 2Vv, where v is the mean atomic 
volume. The value of v may be readily calculated from the composition, 
cell size and density (3-61--0-02 g em~3, determined by flotation) of the 
sample, and is found to be 15-0 4%, The stabilizing ratio would there- 
fore be expected to be slightly less than 1-93 electrons/atom. 
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We have not been able to make the corresponding X-ray measurements 
on a sample of c-(AlFeSi), but can nevertheless calculate the form and 
volume of the Brillouin Zones associated with this compound from the 
lattice parameter (a)=12-523 kx), density (3-61 gcm-3) and composition 
given by Phragmén (1950), together with the knowledge that the phase 
is isomorphous with «(AlMnSi). For if the two phases are structurally 
identical or very closely similar, as indeed they must be to account for the 
observed range of solid solutions reported by Phillips and Varley (1943), 
then their very intense x-ray reflections, and hence the forms of their 
Brillouin Zones, must be identical. We have examined crystals of the 
sample (3/5) of «(AIMnSi) (see Pratt and Raynor 1951 a), and find them 


The first prominent Brillouin Zone of a(AlMnSi). 


to be cubic (a)=12-670+0-000, A, based on Cu Ka, radiation, A= 
1-54050 4) with 138 atoms per unit cell; the atomic arrangement thas 
been determined at least to a first approximation, and will be described 
elsewhere. The x-ray reflections from the {530}, {600} and {352} 
planes are six times stronger than those from any other plane with 
1/2d(hkl) <0-39 A}. a 
The only Brillouin Zone likely to affect the stability of these phases is 
therefore that shown in fig. 2, which is formed by the forty-two {530}, 
{600} and {352} faces. For c-(AlFeSi) these faces lie at p=0-231, 0-237 
and 0-244 4-1 respectively, from the origin of the zone. The volume ce 
the zone (V) is then 0-0612 4-* and the mean atomic volume (v) is 14-9 A®. 


1396 Keith Robinson and P. J. Black on an X-Ray 


The structure of c-(AlFeSi) would thus be expected to be stabilized at 
a composition corresponding to slightly less than 1:83 electrons/atom, 
ie. a slightly lower electron/atom ratio than for the hexagonal 
a(AlFeSi) phase. It is clear that electronically the cubic and hexagonal 
structures do not differ greatly, but other things being equal we should 
expect the cubic phase to have the lower silicon content, since Si has the 
highest ‘ valency ’ of the constituents. 

We may now consider Pratt and Raynor’s suggestion (1951 b) that the 
phases are characterized by a ratio of about two electrons per atom. 
In computing electron/atom ratios from chemical compositions they 
have used the model of Raynor (1944) for the behaviour of transitional 
metals in an electron-rich environment, i.e., ‘effective valencies ’ 
Al—+3, Si=+4, Fe=—2-66, Mn=—3-66. We might also consider 
other possible transitional metal atom valencies, e.g., 

(i) Fe=—2-0, Mn=—3-0, based on the theory of Mott and Jones 
(1936) ; 


(ii) Fe=Mn= 9, as in the f brass structures (e.g. FeAl, Ekman 1931) ; 
(iii) Fe=Mn=-++5-78, as suggested by Pauling (1938). 


Electron/atom ratios calculated from the compositions already given 
would be as follows : 


Fe ‘ valency ’ c-(AlFeSi) «(AlFeSi) 
— 2-66 2-02 2-03 
—2-0 2-14 2°15 

0 2-51 2-53 
+5:78 3°58 3°63 


Estimated from the 
Brillouin Zone ee ate 

These phases may therefore be regarded as electron compounds gover- 
ned by the form of the most prominent Brillouin Zone if each iron atom 
in some way prevents about three otherwise ‘free’ electrons from 
behaving as if free. The Raynor model (Fe valency=—2:66) might 
be satisfactory if the faces closest to the origin of the zones are overlapped. 
Thus a c-(AlFeSi) zone formed solely by the twenty-four {352} faces 
has a volume, V=0-0721 A~%, corresponding to 2-15 electrons/atom. 
For «(AlFeSi), however, the zone formed by {303-10}, {3257} and {3360} 
only raises the possible ratio to 2-01 electrons/atom. 

It is worth noting that the zone calculated from our sample of «(AlMnSi) 
corresponds to an electron/atom ratio about 3°/, lower than that of the 
corresponding c-(AlFeSi), whilst the ratio computed from the composition 
of the sample on the Raynor model is 1-97 electrons/atom, also about 
3% lower than that for c-(AlFeSi). This means that as far as the cubic 
phases are concerned the known compositions and structure are in 
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accord with the view that the compounds are electron compounds in 
which the ‘ effective valency ’ of Mn is algebraically less than that of Fe 
by one unit. 

Moreover, the forms of the most prominent Brillouin Zones of the 
compounds discussed in this paper are consistent with their chemical 
composition only if the transitional metals are given ‘ effective valencies ’, 
if anything, slightly more negative than those suggested by Raynor. 
However, this interpretation assumes that the zones are more full than 
has been necessary in other corresponding cases, e.g. B(AlIMnSi) and 
a(AlCrSi) (Robinson 1952, 1953). If the zones were less full, as in the 
other cases, then the ‘ effective valency ’ of the transitional metal atoms 
would have to be regarded as being even more negative. 


We should like to thank Dr. J. N. Pratt and Professor G. V. Raynor 
for the provision of the specimens examined, and Professor Sir Lawrence 
Bragg, F.R.S., and Dr. W. H. Taylor for their interest in the progress of 
the work. For financial assistance one of us (P. J. B.) wishes to thank the 
Dept. of Scientific and Industrial Research for a Maintenance Grant and 
the other (K. R.) thanks the Ministry of Supply (Basic Properties of 
Metals Committee) and the Royal Commissioners for the Exhibition of 
1851. 
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ABSTRACT 


Measurements of the spectrum of 2-5 Mev neutrons scattered through 
97° from Al, Ni, Bi, W, and Au are presented, together with the 
differential cross section for the processes. 


§1. INTRODUCTION 


Tuts paper is an account of further measurements on the inelastic 
scattering of fast neutrons, using the method already described by the 
author (Poole 1952). The spectrum of scattered neutrons from nickel, 
tungsten, bismuth and gold has been observed, and the measurements 
on aluminium repeated and a more trustworthy result obtained than 
that previously published. 


§2, EXPERIMENTAL METHOD 


For a description of the apparatus, the reader is referred to the earlier 
paper (loc. cit.). 2-5 Mev neutrons emitted at 90° from the D—D reaction 
are scattered by a ring of material placed round the target of the machine, 
and neutrons scattered through (usually) 97° are detected by a 
scintillation counter using an anthracene crystal. A ‘ background’ run 
is made to determine the effect produced by y rays and by neutrons 
scattered from the walls of the room. From the spectrum of pulses 
produced by this counter, the proton recoil spectrum, and so the 
neutron spectrum, may be deduced. The energy resolution is about 
+200 kev at 2-5 Mev. 

Kinematics 

To interpret the curves of scattered intensity the effect of the recoiling 
nucleus must be considered. For scattering through 90° it can be shown 
that the energy of the scattered neutron is 

B,=Hy(M—m)|(M-+m)—E{M|(M-+m)} 
where Hy and £ are the initial neutron energy and energy of excitation of 
the residual nucleus respectively, while M and m are the masses of 
scattering nucleus and neutron. 

Thus if an excitation energy is measured relative to the elastically- 
scattered group (H=0) the energy actually observed will be H{M/(M-+-m)}. 
It is apparent that this correction to HZ is important only in the case of 


the lightest nuclei. For scattering angles less than 90° the effect is less. 
Be Eg BE SP ae 
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§3. RESULTS 
1. Aluminium (fig. 1) 

Groups of scattered neutrons appear at 2-3-L0-15 Mev and 1-4-L0-15 Mey. 
The former is ascribed to elastic scattering of 2-5 Mev neutrons 
(H,=2-32 mev) and the latter to an inelastic group with mean excitation 
of the residual nucleus of 0-9 mev (#,—1-44 Mev). This is consistent 
with the inelastic proton scattering results (Reilley et al. 1952) giving 
levels in *7Al at 0-81 and 0-99 mev, and with the y-spectrum work 


(Hedgran, Hole and Benes 1948, Itoh 1941) giving levels at 0-84 mev 
and 1-05 Mev. 


: Fig. 1 
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Fig. 2 
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2. Nickel (fig. 2) 
Two runs were made with a nickel scatterer. The first showed a group 


with 1-4 mev energy loss, and to make this more clear the experiment 
was repeated with the scatterer placed to receive neutrons emitted at 
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72° to the deuteron beam (H)—2'8 Mev). Although these neutrons are 
no longer strictly monokinetic this shows scattered groups at 2-8 Mev 
and 1-4 Mev, corresponding to elastic scattering, and inelastic scattering 
with 1-4 mev energy loss. This may be compared with inelastic proton 
scattering results of Reilley et al. (1952) giving levels at 1-344 Mev and 
1-48 Mev, and Cowie et al. (1952) who find a level at 1-44 Mev. 
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Spectrum of neutrons scattered from tungsten. 
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Spectrum of neutrons scattered from bismuth. 


3. Tungsten (fig. 3) 

Here apart from the elastically-scattered group at 2-5 Mev, the 
scattered spectrum is continuous, showing a steady rise at low energies. 
Due to poor resolution it is not possible to say exactly what inelastic 
scattering occurs in the region 2-0 Mey to 2-5 Mev scattered neutron energy, 
and therefore the elastic scattering cross section quoted below may 
include a contribution from inelastic scattering with small energy loss. 
This type of spectrum is not unexpected with heavy nuclei, due to the 
probability of a large number of closely-spaced levels. 
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Spectrum of neutrons scattered from gold. 


Table 1. Inelastic and Elastic Neutron Scattering Cross Sections 


Diff. 4n x diff. |Total 
Scatter- cross cross abs. 
Nucleus| E, ing #, section section cross Remarks 
(Mev)| angle (Mev) (barns/ (barns) {section 
sterad.) (barns) 
2:3 0-08 -+0-02 Elastic peak 
27A] 2-5 97° 3:0 
1-4 0:04 +0:02]0-5 +0:25 Probably two close 
levels 
2:8 0-07 +0-02 Elastic peak 
Ni 2°8 70° |_—--— ?__]}——_] 6:0 [——_____ 
1-4 0:05 +0:02/0°6 -+0:25 Probably two close 
levels 
2-5 0-13 +0-04 Elastic peak (plus 
low lying levels) 
W 2-5 | a 720 |———--—-__ 
1:0 to 2:0] 0-10 --0-04]1-3 0-5 Integrated continu- 
ous distribution 
2-4 0:23 --0-04 Elastic peak 
209Bi | 2-4 97° 1:5 0:06 -L0-02 | 0-75-+.0-°3 6:4 | First excited state 
1-1 0:025+0-01]0-3 +0:15 Second excited state 
2:5 0-2 +0:05 Elastic peak (plus 
low lying levels) 
Au 2-5 97° OO 
1-0 to 2:0] 0-07 --0:03]0-9 +0-4 Integrated continu- 


ous distribution 
———E————ES——— | 
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4, Bismuth (fig. 4) 

Groups appear at 2-4+0-1 Mev, 1-5+0-1 Mev and 1-1+0-1 Mev, the 
minima between these being quite low. These may be interpreted as 
elastic scattering together with inelastic scattering from excited states at 
0-9-L0-2 mev and 1:3-40:2 mev. The level spacing is unusually high for 
a heavy nucleus, as would be expected from the nearness of *°°Bi to the 
‘magic’ numbers for both neutrons and protons. 


5. Gold (fig. 5) 
Shows elastic peak at 2-5 Mev plus a continuous distribution of neutrons, 
rising at low energies, consistent with many closely-spaced levels. 


6. Cross section values 


Values for ‘light’ elements (discrete levels) are differential cross 
sections for scattering through the stated angle, leaving the nucleus in 
the level concerned. Values for ‘ heavy ’ elements (continuous spectrum) 
are split into above 2-0 mev (‘elastic ’ scattering) and 1-0 Mev to 2-0 Mev 
(‘inelastic ’ scattering). Additionally for comparison the integrated 
cross section assuming a spherically symmetrical angular distribution, 
and the total absorption cross section (taken from the literature) have 
been tabulated. Spherical symmetry is probably only justifiable in the 
case of inelastic scattering from many levels (‘ heavy’ nuclei). 
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ABSTRACT 


The optical method of Burge and Jackson (1951) is employed to deter- 
mine the rate of extension of the helium II film over a previously bare 
solid surface. The passage of the front of the film over two regions of 
the surface a known distance apart is observed, and a provisional 
apparatus has shown that the method is feasible and has indicated what 
changes should be made to obtain significant results. 


§ 1. INTRODUCTION 


FRENKEL (1940) has derived a formula for the rate of the initial spread 
of a film of a normal liquid over a vertical solid surface partly immersed 
in the liquid. For the case of water the velocity of the outer layer of the 
film at a height of 1 cm above the liquid surface is calculated to be about 
2mm/s. In the case of liquid helium just above the A-point, the corres- 
ponding value will be about 80 cm/s. The formula cannot, of course, 
be applied to the case of helium II because of the superfluid properties 
of the latter. No attempt has yet been made to derive a formula appro- 
priate to the latter case. It therefore seems of interest to make a direct 
experimental determination of the rate of extension of the helium II 
film over a previously bare solid surface. 

It seemed probable, from the observations of Burge and Jackson (1951) 
of the helium IT film by an optical method, that it should be possible to 
observe the passage of the front of the film over two regions of the solid 
surface a known distance apart and so derive the velocity of extension of 
the film. Accordingly a provisional apparatus was devized, which 
involved the least possible alteration to the present optical system, and 
constructed from available material. This apparatus has shown that the 
method is feasible and has indicated what changes should be made in a 
future version to obtain accurate and reproducible results. As we do 
not anticipate being able to carry out observations with the new version 
of the apparatus in the immediate future, it is thought to be of interest 
to publish this preliminary account of the work. 


* Communicated by the Authors. 
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§2. APPARATUS AND MuTHoD 


The apparatus is shown in fig. 1. It consists essentially of two small 
stainless steel mirrors, M, and Mg, each 1 cm long, one above the other in 
the same vertical plane and connected by a wire of diameter 1mm and about 
30 cm long coiled into a spiral. As actually constructed the two halves 
of the coil were wound in opposite directions to avoid the mirrors being 
displaced from the same plane as the apparatus cools from room tempera- 
ture. The arrangement is supported by a bent nickel rod N, the upper 
end of which is connected to a thin-walled glass tube T terminating in 


Fig. 1 


the pin-chuck shown above and the lower end is attached below the lower 
mirror at B. As in the experiments of Burge and Jackson, half of each 
mirror is coated with a layer of barium stearate one moleeule thick and 
the other half with a layer three molecules thick. Sodium light linearl 

polarized at 45° to the plane of incidence is reflected from both ore 
and is observed through a 4/4 plate and a nicol prism. If the two latter 
are so adjusted that the intensity of the light reflected from the two halves 
of the mirrors is the same in the absence of the helium II film, then th 

passage of the film over a mirror can be observed by the lack of meant 
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of the ‘1’ and ‘3’ areas of the mirror. The experiment consisted in 
principle in dipping the tip A into liquid helium IT and measuring the 
time interval between the appearance of the film on M, and M,. 

Dipping the apparatus directly would have disturbed the optical 
settings and so a cylindrical spring-loaded plunger P was arranged with 
a central hole in it, operated externally, so that having set A 1 mm or 
so above the liquid level, contact could be made by lowering the plunger 
until it was immersed and so raising the level of the liquid above the 
tip A. The whole apparatus was surrounded by a cylindrical copper 
radiation shield 8, which had appropriate windows cut in it and was 
supported by another glass tube which joined T at C just below the 
pin-chuck and about 20 cm above the mirrors. Since the chuck was in 
good thermal contact with a metal vessel containing liquid hydrogen at 
20°K, the temperature of C would always be above the A-point, and so 
the mirrors would not be covered with the film by the latter creeping 
over the surface of S, up the glass support and down T and N. 

To obtain significant results it was essential that there should be no 
appreciable temperature gradient in the mirror system. Enclosing the 
latter inside the well-conducting copper shield, which was always in 
contact with the liquid helium and was covered by the helium II film, 
largely satisfied this requirement, but the method of supporting the 
mirror system was found to result in too large a heat influx (see later), 
and this must be improved in a future form of the apparatus. 


§ 3. RESULTS 


It was essential for the purpose of the experiment that a film should 
not form from the vapour phase to a thickness comparable with that 
from the liquid phase during the time of observation. This was first 
checked in a preliminary experiment. 

Observations were first made immediately below the A-point but no. 
films were detected on M, or M, until the temperature was reduced to. 
2-09°x when A was in contact with the liquid. It then appeared that the 
heat influx down the wire support was sufficient to evaporate the film on 
A before it reached the junction at B for any temperature between 2-09°K. 
and 2:19°x. For any temperature below 2-09°K the creep rate was suffi- 
ciently great for a film to pass B in spite of the heat influx, which was. 
estimated to be about 100 erg/sec. A reduction of this figure by a change 
in the method of supporting the mirror system is clearly desirable in any 
future apparatus. 

Observations were next made at 1-90°x and the passage of a film over 
M, and M, could be observed. It was, however, found necessary to fix 
the attention on one mirror for some 10-15 sec before it was possible to 
detect the subsequent lack of matching due to the passage of a film over 
the mirror. It was therefore not possible to observe the movement of 
the film past M,, transfer the attention of the observer to M, and detect. 
the movement of the film over the latter if the interval was less than 
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15 sec. However, by concentrating attention on M, the following 
intervals were observed between the time of contact of the liquid with A 
and the arrival of the film on M,: 4:5, 4:2, 4:5, 4-5, 4:5, 4:0, 4-7 sec, with 
an average of 4-4 sec. 

The film did not appear as a sharp edge passing over the mirror, but as 
a gradual change in the relative intensities on the two sides of the 1-3 
boundary lasting for less than 1 sec which is included in the figures given. 
This was to be expected since the film was flowing on to the mirror from 
the much smaller periphery of the wire and having to ‘ build up ° to its 
equilibrium thickness. 

A repetition of the observations at 1-44°K gave the following times : 
3:5, 4:3, 4:3, 3-5, 4:3, 4:5, 5-3, 3:5, 4:3, 5-5, 4-4, 4:0, 44sec; average 
4:3 sec. Each time contact with the liquid was broken the film evapor- 
ated from M, in about 5 sec, or rather a considerable reduction in thickness 
of the film took place as no attempt was made to measure the actual 
thickness in these experiments. 

It would therefore appear that when further developed the apparatus 
will be capable of giving an accurate figure for the rate of extension of 
the helium II film over a solid surface. No actual velocity will be calcu- 
lated from these preliminary observations because of the uncertainties 
due to (1) the finite time of spread of the film from the wire to the larger 
periphery of the mirrors, and (2) the fact that the film divides at B. It 
is not known whether the film preserves the same thickness on the mirror 
system and the corresponding point on the wire support N and so travels 
at a slower rate than on the wire at A, or whether the velocity remains 
constant and the film thickness decreases. 

In the new apparatus now being constructed the method of support 
is being changed so that the heat influx is reduced and the division of the 
film at B is avoided. By making the mirror system of a continuous 
flat polished strip the change of periphery at the mirrors is avoided, and 
the latter are arranged side by side in the same vertical plane so that the 
two 1-3 boundaries are co-linear and no movement of the observer’s 
eye from one mirror to the other will be required. 
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CXLIX. CORRESPONDENCE 


Estimation of Mean Lifetimes from Multiplate Cloud Chamber Tracks 


By M. S. Bartnerr 
University of Manchester* 


[Received October 14, 1953] 


In a recent discussion (Bartlett 1953) of the statistical estimation of the 
mean lifetimes of unstable particles, the methods proposed were illustrated 
with reference to some data of Fretter et al. (1953) on neutral V particles. 
This may prove confusing to the reader, as I had overlooked the fact 
that these authors were making use of data from a multiplate cloud 
chamber, and the estimation equations are not strictly applicable in this 
case without modification. My attention was kindly drawn to this 
point by J. A. Newth (cf. Page and Newth 1954), and while the appro- 
priate modification of the estimation equation has already been indicated 
and used by Bridge et al. (1953, § TX), it may be helpful to record it here, 
together with the further modified formulae if the accuracy of the estimate 
is to be precisely evaluated. 

In relation to the track of the rth of n particles observed to decay in 
the chamber, there will now be, in addition to the time of decay ¢,, measured 
from an appropriate zero and the total possible time 7,, a number of 
further times 7',, (s=1, 2,...) when the track or its continuation enters 
and leaves the intermediate plates (if any). It is convenient to write 
as before ¢,/0=u,, T',/0=U,, where @ is the unknown mean lifetime ; we 
now have also 7',,/9=U,,._ I shall use the further notation 

A,=1—e74%+2,(—1)* e~ Us, 
B= ee oF fe ae Ni) Able e- ora, 
C,=U,2 0% —3,(—1)U 2 0 M8 5 
faa iA Welles 
Pad ia a Lede 
d,= Ta | ain Lee. 
Then the probability distribution of the decay times ¢, is 
n n aU. 
I] er dt,= I ee 51) 
r=1 r=1 r 
and the estimation equation (2, denoting summation over the n different 
tracks) is 


02, 5 =3,(14y—14B,/4,)=0. aera tins. (2) 


* Communicated by the Author. 
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Moreover, 


aL OF. BB 
wns, (— Fe) fain a (I- Fg) = 


the approximate standard error of the estimate of @ obtained from (2) 
being 6/1/2,(I,),, or equivalently that of 1/@ being 1 |[0/%,(L),]- 

The standardized quantity S(9) with zero mean and unit standard 
deviation, most suitable for plotting (against, say, convenient values of 
1/0) if eqn. (2) is to be solved by interpolation or confidence limits for @ 
are required, is 


SO)=2,(4,—1+ BJA) 20) 

The higher moments of S(#) could be worked out if required, though they 

will now be rather complicated. Even for the evaluation of 2,(/,), and 

S(0) we cannot unfortunately make much use of the table provided in 

my paper; direct evaluation for one or two simple values of 1/6, while 
tedious, is, however, still quite practicable. 

In the particular region 1/9 ~ 0, the expression for S(#) becomes 

X,(t,—3C,/6,) 1 

S(@) ~ ine EX eae aes . . . . . (5) 


where d,, C, \2 
x= x— (i) 


Thus for 1/6 ~ 0, the estimate becomes 

1 X(t, — $¢,/0,) . 

eb with standard error 1/1/2,X,, EY 
(the negative sign outside the expression on the right was inadvertently 
left out of the corresponding eqn. (6) of my paper). This rough 
solution must be used with caution, for it depends on 7',/0 being small; 
it is not recommended if S( co) is significantly negative, say less than —2. 
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Crystallographic Slip in Diamond 


By M. Sean and J. W. MentErR 


Research Laboratory for the Physics and Chemistry of Surfaces, Department 
of Physical Chemistry, University of Cambridge 


[Received October 1, 1953] 


THE technique of reflection electron microscopy developed by von 
Borries ( 1940) and more ‘recently by Menter (1952, 1953) and others has 
been applied to an investigation of surface features on polished diamonds. 
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In this method the surface is illuminated with electrons incident at a 
glancing angle of 1-2° and viewed at an angle of 6-12°. Contrast arises 
primarily from shadows cast by asperities standing out from the surface. 
Thus the most prominent features seen on polished diamonds examined 
in this way are parallel grooves formed by the abrasive. (All the diamond 
surfaces examined had been polished in one direction only.) The 
visibility of the grooves is greatest when they are oriented at right angles 
to the plane of incidence but surface features within the grooves are not 
seen on account of the shadows cast by their edges. However, the 
grooves may be rendered almost invisible by arranging to have them - 
parallel to the plane of incidence. Detail within the grooves then 
becomes more readily visible and it has been found that there are a number 
of roughly parallel lines on the surface, different in character and direction 
from the polishing marks. These lines can be seen more easily by 
optical examination of the carbonaceous contamination layer which is 
formed on the surface by the action of the electron beam. The nature 
and origin of this contamination has been investigated by Ennos (1953). 
He showed that it is formed only where the electrons strike the surface. 
If, therefore, a specimen is irradiated in the electron beam and sub- 
sequently examined in the optical microscope, it has the appearance of 
having been shadowed at a very oblique angle. 

Figure 1 (Pl. 52, figs. 1 and 2) shows an optical micrograph of a polished 
diamond surface which has been contaminated in the electron microscope. 
The directions of the polishing marks and of the incident electron beam are 
indicated by the arrow (A). There are also several cross lines (B) 
rendered visible by the contamination shadowing effect. Such lines 
correspond to changes in surface level which may in some cases be as 
small as 20 A. This has been established by stripping the contamination 
layer from the surface and examining it at high magnification as a low 
angle shadowed replica in the transmission electron microscope. The 
most likely explanation of these lines seems to be that they are due to 
crystallographic slip occurring during polishing. Extremely high local 
pressures must be produced whilst polishing and these may well be large 
enough to initiate plastic deformation (the Knoop indentation hardness 
of diamond is 5 500-6 950 kg/sq. mm). 

Three polished diamonds have been examined. They are all of type I 
and x-ray and electron diffraction photographs show no evidence of 
macroscopic twinning. In each case the mean direction of the cross 
lines follows the trace of a {111} plane. There is, however, a spread of 
several degrees in their direction and this could be due to cross-slip. 
The lines only follow one set of {111} planes on each stone and, asssuming 
a glide direction [110], these are the octahedral planes of greatest resolved 
shear stress (Schmid criterion). On one diamond the face which had been 
examined was repolished in the same direction as before, the thickness 
of the layer removed being less than 25. Hxamination of the same 
field before and after this second polish showed that, although the lines 
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were visible in each case, there was no detailed correspondence between 
them. The stone was then repolished again but this time in a completely 
different direction. Lines were visible on a new set of {111} planes and 
the old cross lines could not be detected. It would seem unlikely, there- 
fore, that the lines could arise from some lamellar growth structure within 
the diamond or from cracking of the surface. Moreover, the lines do not 
follow the traces of octahedron planes of greatest resolved normal 
breaking stress (Sohncke criterion). 

Further confirmation that the lines are not due to cracking has been 
obtained by examining a gross scratch produced on a polished dodeca- 
hedron face by a hard octahedron apex (fig 2). The gross white lines 
extending a short distance from the scratch follow the {111} planes of 
greatest breaking stress and are certainly cracks. The second set of lines 
(B) are the polish marks and there is a third set (C) following traces of 
another set of {111} planes, which could be slip lines. 

Williams (1932) has suggested that some of the lines observed on 
natural diamonds might be due to slip and Tolansky and Omar (1953) 
have described a case of possible slip on opposite octahedron faces of a 
twinned diamond, as mined. The present observations suggest that 
plastic deformation may be produced in the act of polishing diamond. 
This would have an important bearing on theories of the mechanism of 
the polishing of diamond. 
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Optical micrograph of contaminated diamond surface. A—=direction of polishing 
marks and incident beam, B=lines probably due to slip. 850 


Fig. 2 


Reflection electron micrograph of gross scratch on diamond, B=polish marks, 
C=lines probably due to slip. 
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CL. Notices of New Books and Periodicals received 


The Stability of Rotating Liquid Masses. By R. A. Lyrtueron. [Pp. vii+-150.] 

(Cambridge University Press.) Price 35s. 
THE forms of gravitating fluid masses in steady rotation and their stability 
have been the subject of many monumental papers. The present book explains 
the leading principles and results, while referring to the original papers for some of 
the details of the calculation. The theory of the exchange of stabilities, developed 
by Poincaré for these problems, has much wider application, and the account 
of it here is the clearest yet given. The account of ellipsoidal harmonics, 
though much shorter than those in the original papers and in Hobson’s book, 
is full enough for the purpose and should be useful to specialists in other fields. 

The outstanding results of the theory are that of Poincaré, that the triaxial 
ellipsoid becomes secularly unstable at a certain speed of rotation; that of 
Liapounoff and Jeans, that when this speed is slightly exceeded there is no 
neighbouring configuration that is secularly stable ; and that of Cartan, that 
when the Jacobi figure becomes secularly unstable it simultaneously becomes 
ordinarily unstable. Jeans had overlooked this possibility, saying that secular 
instability must always occur first. The importance of Cartan’s result was 
first noticed by Lyttleton. Ordinary instability implies that a small dis- 
turbance will increase rapidly without oscillation, and further examination 
shows that the only possible final state is one where the fragments separate 
with a hyperbolic velocity. A star acquiring too high a rotation as it condenses 
cannot become a binary star; it would become two independent stars. The 
fission theory of the formation of binary stars must be abandoned. 

The original papers, apart from their unavoidable difficulty, are obscure in 
some places and inconsistent in others, and Lyttleton’s careful study of their 
interpretation is a valuable clarification. Fad; 


Luminescence and the Scintillation Counter. By S. C. Curran, F.R.S. 
(Butterworths Scientific Publications, London, 1953.) [Pp. 205-+xciii.] 
Price 32s. 6d. 

THE extensive development and use of the scintillation counter in nuclear 

physics has for some time warranted a detailed treatment of the counter and 

its principles. The present volume amply fulfills that need. The author is well 
qualified for the task. He was one of the first physicists to ‘ modernize the 
spinthariscope by inclusion of a photoelectric multiplier to record the emission 
due to single alpha particles. One of the features of the text is that it treats 
both the operational characteristics and the fundamental processes in each 
component of the counter in some detail, prefaced by a valuable chapter on the 
interaction of the various nuclear particles and high energy radiation with 
matter. Such subjects as limiting sensitivity and time and energy resolution 
of the counter and the detailed characteristics of inorganic and organic phos- 
phors, inclusive of liquid and plastic scintillators, are given full discussion. 
A comprehensive though brief account of applications of the counter to radiation 
detection and energy measurement is given in a later chapter. As additional 
topics the text includes a description of the application of photomultipliers to 
direct detection of particles and a treatment of some electrical circuits involved 
in the use of scintillation counters, such as those for short time resolution, 
fast scaling and pulse magnitude discriminators and analysers. Full and up 
to date references are given at the end of each chapter and name and subject 
indices at the end of the text. There is no doubt that this volume will be 

a valuable asset to every nuclear physics laboratory, since within its compass it 

includes information which up to the present has been widely scattered in 

original papers. Its value is also enhanced by the proper connection it makes 
between the physics of constituent parts of the scintillation counter and the 
application of the whole, a conjunction which will greatly assist in assessing the 

possibilities of the counter for any specific use. G. F. J..G. 
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Modern Radiochemical Practice. By G. B. Coox and J.F.Dunean. (Oxford: 

Clarendon Press, 1952.) [Pp. xx+407.] Price 42s. 

Tuts book by two members on the staff of the Atomic Energy Research 
Establishment at Harwell is a good example of the impossibility of separating, 
in the science of radioactivity, the spheres of interest of chemists and physicists. 
It is called Modern Radiochemical Practice, but a very considerable part is 
devoted to the discussion of subjects which a chemist of the old school would 
consider as pure physics, as, for instance, the radiations emitted by radioactive 
material, the methods of their measurement and the determination of their 
energy, and the description of very special electronic equipment. However, 
the authors are certainly right in assuming that no laboratory interested in 
radiochemical work can be confident of obtaining reliable results unless at least 
one of the workers—perhaps, best of all, a physicist—is fully acquainted with 
the theory and practice of the measuring apparatus, and nobody will, therefore, 
object to the inclusion of this material, even if the title of the book seems to 
restrict its contents to ‘radiochemical’ practice. (The further qualification, 
‘modern ’, might well have been omitted ; surely no one would suspect that a 
publication by two Harwell men and the Clarendon Press, appearing in first 
edition in 1952, might be old-fashioned.) 

The book is meant principally as a guide for those who want to apply 
‘ radioelements as indicators ’, as one used to say in earlier days, or ‘ isotopes 
as tracers’, as the modern parlance prefers. The substitution of the word 
‘tracer’ for ‘indicator’ seems to be an advantage because the new term is 
more restricted and precise ; on the other hand the meaning of ‘ isotopes ’ as 
understood by the authors is hardly to be recommended as it is bound to lead 
to contradictions. This technical term was coined by F. Soddy to denote two 
different species of atoms belonging to the same element and thus occupying 
‘the same place’ in the Periodic System ; but the authors have adopted the 
recent loose habit of calling every radioactive species an isotope ; so for instance 
when they append to the text of their book a table, not of ‘ atomic species ’ or 
‘nuclides’, but of ‘isotopes’, or when they speak of ‘isotopes of the same 
element ’. In other places, however, they are quite aware of the strict meaning 
of isotopes ; e.g. if they discuss the ‘ separation of isotopes’, or explain the 
method of ‘ isotopic dilution ’. Let us hope that the term ‘ nuclide ’, now widely 
accepted in the U.S.A., will also in the English literature replace the use of 
‘isotope’ in the broader sense, thus restoring this indispensable term to its 
accurate original meaning. 

A special feature of the book is the last chapter, ‘ Suggested experiments with 
radioactive material’. Here thirty very different experiments are described 
in sufficient detail for the introduction of beginners who will be grateful for the 
many practical hints included. It is obvious that the wide experience of the 
authors has enabled them to produce a book which will be found very valuable 
by scientists in many fields who intend to use radioelements as tracers. 
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Fluorescence of Solutions. By E. J. Bowen and F. Woxgs. (Longmans, 
Green & Co., Ltd.) [Pp. 90.] Price 25s. net. 


Tats book takes the form of a useful monograph for students and practical 
workers who need an elementary introduction to the subject. It avoids 


detailed mathematical proofs and devotes considerable space to the methods 
and practical side of measurements. 
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